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Abstract 

We derive new results about properties of the Witten zeta function 
associated with the group SU( 3), and use them to prove an asymp¬ 
totic formula for the number of n-dimensional representations of SU( 3) 
counted up to equivalence. Our analysis also relates the Witten zeta 
function of SU (3) to a summation identity for Bernoulli numbers dis¬ 
covered in 2008 by Agoh and Dilcher. We give a new proof of that 
identity and show that it is a special case of a stronger identity involv¬ 
ing the Eisenstein series. 


1 Introduction 

1.1 Asymptotic enumeration of representations 

Let p(n) denote the number of ?i-dimensional representations of the group 
SU( 2), counted up to equivalence. Since SU( 2) has (up to equivalence) 
one irreducible representation 14 of each dimension k = 1 , 2 ,..., each n- 
dimensional representation ©£1 is encoded by a unique integer parti¬ 

tion n = Ai + ... + A m , such that rp is the number of parts in the partition 
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that are equal to k. Thus pin) is the number of integer partitions of n. For 
large n, p{n) is given approximately by the well-known Hardy-Ramanujan 
asymptotic formula [20J (see also [22, Chapter 2]) 

p(n) = (1 + o(l))— \—E K '/ 2n ! 3 as n -» oo. (1) 

4-\/3n 

This work began with the idea of answering a similar question for the 
group SU{2>). Let r(n) denote the number of n-dimensional representa¬ 
tions of the group SU(3), counted up to equivalence, or equivalently the 
number of n-dimensional complex-linear representations of the Lie algebra 
A ‘2 = sl{ 3, C). The irreducible representations of 577(3) are a family of 
representations Wj.k indexed by pairs of integers j,k > 1, where it is well 
known (see, e.g., Cl Chapter 5]) that dim Wj^ = \jk(j + k ). A general n- 
dimensional representation decomposes as a sum of these Wj^ s, each with 
some multiplicity. Thus, it is easy to see that the numbers r(n) satisfy the 
generating function identity 


E r (") x "= II 

n =0 j,k=1 


( 2 ) 


(with the natural convention that r(0) = 1) analogous to Euler’s product 
formula £“ =0 p(n)z n = YYZLii 1 — x k ) 1 for the generating function of parti¬ 
tions. Using ([2]) one may easily compute the first few terms of the sequence 
r(n), which are (starting at n = 0) 1,1,1, 3, 3, 3, 8, 8, 9, 17, 19, 21, 35,39,.... 

One of our main results is the following analogue for the group SU (3) of 
the Hardy-Ramanujan asymptotic formula ([Tj). 

Theorem 1.1. As n —>• 00, we have 

r{n) = (1 + o(l)) exp (d^ 2 / 5 - A 2 n 3 / 10 - A 3 n 1/5 - ^n 1 / 10 ) , (3) 

for constants K, A \, A 2 , A 3 , A 4 , which are defined as follows. LetT(s) and 
£(s) denote the Euler gamma function and the Riemann zeta function, re¬ 
spectively. Denote 


x = (l r (i) 2 C (t)) 3/U \ 
y =-V5FC c (i). 


(4) 

(5) 
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Then A\, A 2 , A%, are given by 


Ai 

= 5X 2 

= 6.858260476163126..., 

(6) 

A2 

II 

* 

1 

= 5.773601745105114..., 

(7) 

^3 

II 

0 w 
* 

to 

= 0.911341072572436..., 

(8) 

a 4 

= 11 x~ 7 y 3 

3200 

= 0.351637541558209..., 

(9) 


and the multiplicative constant K is given by 


K= ^ x ^ 

V5 


X~ 10 Y 4 ) = 2.4462903348641789 ... (10) 


2560 


The formula (10) for K is a simplified version of a more complicated for¬ 


mula that appeared in the preprint version of this paper. The simplification 
depends on a result from a recent paper by Borwein and Dilcher |12| —see 


Theorems 1.4 and 1.5 and the remark following Theorem 1.5 below. 


1.2 The zeta function of SU(3 ) 

Our efforts to derive and prove the asymptotic formula ([3]) led us to make 
several additional discoveries that are of independent interest. Understand¬ 
ing the asymptotic behavior of the sequence r(n) turns out to require a 
close study of a natural Dirichlet series (or zeta function) associated with 
the representation theory of SU(3), namely 

OO OO . 

<n) 

This function has an interesting history dating back to the 1950’s papers 
of Mordell [33] and Tornheim [45]. It is (up to the trivial factor 2 -s ) a 
special case of a so-called Witten zeta function, which is more generally 
a Dirichlet series C 0 (s) associated with a semisimple Lie algebra g and defined 
by 

CoM = 

p 

where the sum extends over the irreducible representations of g. (The same 
function in the context of groups has been referred to by some authors as 

the representation zeta function; see, e.g., [3], pi].) 
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Much of the work to date on u>(s) has concerned its evaluation at positive 
integers. This was motivated to a large part by Mordell’s discovery [33] that 
the even integer values uj(2n) are rational multiples of 7r 6n , being given 
explicitly by the formula 

(—1 ) n + 1 f97r) 6n f 1 

w(2n) = — ' } / B 2 n (x) 3 dx (n > 1), (12) 

6((2 n)\y J 0 

where B m (x ) denotes the mth Bernoulli polynomial. Later authors have 
since extended Mordell’s results in several ways. Witten [33] related the 
values at even integers £ g ( 2 n) (n > 1 ) of the zeta function of a general 
semisimple Lie algebra 0 to the volumes of certain moduli spaces of vector 
bundles of curves. His formulas imply that £ g (2n) is a rational multiple of 
7 r 2rn , where r is the rank of the Lie algebra g. Subbarao and Sitaramachan- 
drarao [33] proved another explicit formula for u>(2n), namely 

w(2n) = -J^f JC(2fc)C(6n- 2k) (n > 1). (13) 

According to Zagier’s paper [50j, the same formula was apparently redis¬ 
covered by Garoufalidis and Zagier and independently by Weinstein. In the 
same paper, Zagier also gives without proof an analogous formula for the 
odd integer values u(2n + 1), namely 


u>(2 n + 1) = —4 


k =o 


4n — 2k + 1 
2 n 


C(2fc)C(6n-2fc + 3) (n > 0). 


The same formula as well as its proof was published a short time later by 
Huard, Williams and Nan-Yue [22], who seem unaware of Zagier’s paper. 
The formula (13) is also given in the paper of Gunnels and Sczech [18] along 


with an analogous formula for the values of the zeta function of <SC/(4) at 
even positive integers. 

As an illustration of the above formulas, here are the values oj(m) for 
1 < m < 5: 


w(l) = 2C(3), 


w(3) = -2vr 2 C(7) + 20C(9), 

19 12 

“ 273648375^ ’ 
o 70 

w(5) = -g A(H) - yvr 2 C(13) + 252C(15). 
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Several of the papers mentioned above also contain formulas of a similar 
flavor for the values of more general “multiple zeta values” defined in terms 
of several integer-valued parameters; see also [T7] , [23], [35] . [5T] for related 
results. 

The results described above may not make evident why it is natural to 
consider w(s) as a true Dirichlet series (that is, as a function of a complex 
variable s), but plenty of precedents from the history of analytic number 
theory suggest that this is worth doing. Perhaps guided by such considera¬ 
tions, Matsumoto [25] initiated the study of the analytic continuation of the 
three-variable Dirichlet series 

OO 

(mt(si,s 2 ,s 3 ) = j~ Sl k~ S2 {j + k)~ S3 , (14) 

j,k=l 


which contains w(s) as the specialization si = S2 = S3 = s, and several 
other multiple Dirichlet series of a somewhat similar nature. He referred to 
Cmt(si, S2, S3) as the Mordell-Tornheim zeta function. He and his col¬ 
laborators later extended such results to more general (multivariate) Witten 
zeta functions m, eh, m, 12a. 

Following in the footsteps of Matsumoto’s pioneering work, and guided 
by a specific need to understand the complex-analytic properties of co(s) in 
connection with our asymptotic analysis of the sequence r(n), we proved 
several results about w(s). The first one is as follows. 

Theorem 1.2. (1) The series © converges precisely for complex s£C 
satisfying Re(s ) >2/3, and defines a holomorphic function in that region. 

( 2 ) u(s) can be analytically continued to a meromorphic function on C. 

(3) uj(s) possesses a simple pole at s = | with residue 2 J^ T Q) 3 , and 

for each k = 0, 1 , 2 ,... it has a simple pole at s = \ — k with residue 
(—l) fc 2 _ 4 fc ( 2 fc')C -ft has no other singularity points. 

(4) | w(s)| grows at most at a polynomial rate for fixed Re(s ) as Im(s ) —> ±00. 
More precisely, for any closed interval I = [a, b] C M, there exist constants 
C,M > 0 such that the bound 

\u(a + it)\ < C\t\ AI (15) 


(where i = y/—l) holds for all a € I and |t| > 1. 


We note that Matsumoto already proved part (2) of Theorem 1.2 (see 
Theorem 1 in his paper |25]), but he was studying the more general function 
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(14) of several complex variables and did not examine more closely the 


properties of cu(s). 

At the end of his discussion of the Mordell-Tornheim zeta function, 
Matsumoto writes: “It is now an interesting problem to study the prop¬ 
erties of the values of Qmt{si, $ 2 , S 3 ) at non-positive integers .” We solve 
this problem for the case of our function u(s) = (mt{s, s, s ). The result is 
as follows. 


Theorem 1.3. oj(s) has the values 

cu(0) = 
co(—n) = 0 


(n = 1,2,3,...). 


(16) 

(17) 


Note that the equations ( [l6| ) -© imply that Mordell’s property that 
u}(2n) is a rational multiple of tt 6 " holds true for all integer n. Furthermore, 
somewhat surprisingly, the formula ( |13[ ) of Subbarao and Sitaramachan- 
drarao still holds true for n = 0 if one accepts the convention that (~ j) = 1, 
even though their proof clearly does not apply in that case. (Perhaps less 
strikingly, the formula also holds true for negative values of n, simply because 
the range of the summation becomes empty.) By comparison, Mordell’s for¬ 


mula (12) gives the incorrect value of —1/6 in the case n = 0. 


Our final result on cj(s) is a formula expressing its derivative at s = 0 in 
terms of familiar special constants and a curious definite integral. 

Theorem 1.4. Let 7 denote the Euler-Mascheroni constant. The value 
u/(0) is given by 


1 . 3, _ 1 C (| + tt) ( (—| — it) 


cu'(0) = ^(1 + 7) + l log(27r) - 2C'(-1) + 2 / 
= 1.83787706640934548356.... 


(| + it) cosh( 7 rt) 


dt 


(18) 


Note that the integral in (18) converges very rapidly due to the expo¬ 


nential decay of the factor l/cosh(-7rf), and is easy to evaluate numerically 
to high precision. Curiously, the numerical value of the integral (including 
the multiplicative factor |) is —0.002807659..., so its relative influence on 
u/(0) is very small—less than a fifth of a percent. 


Update added in revision. Following publication of the preprint version 
of this paper, the following much simpler formula for u/(0) was proved by 
Borwein and Dilcher m- 
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Theorem 1.5. u/( 0 ) = log( 27 r). 


See also [9] for related recent developments. 

Theorems |1.2f T.5 will all play a role in our asymptotic analysis leading 
up to the proof of Theorem 0 In particular, the poles at s = 2/3 and 
s = 1/2 and their residues, and the evaluation w(0) = /, have a large 
influence on the form of the asymptotic formula ([3]). The evaluation of 
u/(0) enters into the definition of the multiplicative constant K, which ends 
up emerging from our calculations in the form 


K = 


VI 

VVr 


x 1/3 


exp 


2560 


X- 10 Y 4 + w'(0) 


(19) 


(clearly equivalent to (10) in view of Theorem |1.5|). The poles with negative 


real part and the location of the zeros © do not enter into our present 
analysis of the asymptotic behavior of r(n), but those results, aside from 
being interesting in their own right, will also start having an effect on the 
asymptotics of r(n) if one should attempt to derive more detailed expansions 
for r(n) which include additional lower order terms currently encapsulated in 
the (1 + o(l)) term in ([3]); see the discussion of open problems in Section 12 


1.3 The Bernoulli numbers and Eisenstein series 

One of our surprising discoveries was that Theorem |1.3| has an intriguing 
connection to the Bernoulli numbers, and ultimately to the theory of modu¬ 
lar forms. Recall that the Bernoulli numbers are defined by the power series 
expansion 

r °SL r 

X _ \ n 

e x — 1 ^ n\ 

n =0 

and that they are related to special values of the Riemann zeta function by 
the well-known equations £(1 —n) = (—1 ) n+1 ^f, C{2n) = (—l) n+1 ^ 2 n, 

(n > 1). 

It is well-known that the Bernoulli numbers satisfy a host of remarkable 
summation identities (see 0 , m , and im pp. 81-85] for many examples). 
Our analysis of the properties of ui(s) led us to discover one such identity, 
namely the statement that for n > 1 we have 


Bqu+2 _ _ (4n + 1)! / 2 n \ E>2n+2k _Jhn-2k+2_ 

6n + 2 “ (2n)! 2 ^ V 2 ^ - V 2n + 2k ’ An-2k+ 2' 


( 20 ) 
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Equivalently, in terms of the values at integer arguments of the zeta function, 
(20) can be rewritten in either of the forms 


/-la | r)^ (4n+l)!^ Gfc-l) , r,,\slA Of , o^ 

C(6n + 2) — , 1 • /0 _ u2 7 67 yC(2n + 2fc)c,(4n — 2fc + 2), 


C(-6n- 1) = 


6n+1 (2 „)P 

(4n + 1)! 


( 21 ) 


2 n 


(2 n )! 2 ^ V 2 fc -1 

v ' fc=i 


E 


C(-2n - 2k + l)C(-2n + 2fc - 1). 


( 22 ) 


The identity (20) turns out to have been discovered and proved earlier, be¬ 
ing due to Agoh and Dilcher in their 2008 paper [2]. It belongs to a class 
of several summation identities that Agoh and Dilcher refer to as lacu- 
nary recurrences, since they have the appealing property of expressing a 
given Bernoulli number in terms of a relatively small number of Bernoulli 


numbers of lower index; in the case of (20), Bq u+2 is expressed in terms 
of -£? 272 + 2 f -^ 2 n+ 4 ; ■ • • ? -^ 4 n* See also m for further discussion of lacunary 
Bernoulli number recurrences. 


Note that (20) does not hold for n = 0, when the left-hand side is equal 


to 1/12 and the right-hand side is equal to 0. 

Being at first unaware of Agoh and Dilcher’s result, we independently 


proved (20). This was fortunate, since our method of proof actually implies 


a stronger result that does not follow from the methods of [I], and gives 
hints of a deeper theory at play. To formulate our new result, recall that 
the Eisenstein series of index 2k is an analytic function defined for 2 

in the upper half-plane M = {x + iy : y > 0} by 


G 2k (z ) = 


E 


1 


(m,n)gZ 2 \{( 0,0)} 


(■mz + n) 2k 


(k > 2). 


The function G 2 k(z) is related to the zeta value ('(2k ) (hence also to the 
Bernoulli number B 2k ) by the well-known relation 


lim G 2k (z) = 2C(2 k), 

Im( 2 :)—>-oo 

usually formulated as the statement that 2</(2 k) is the constant term in the 
Fourier series G 2 k{z) = 2((2k) + Y^=i c 2 k(n)& 2mnz (see [7] for these facts 
and additional background on modular forms mentioned below). We are 
now ready to state our result. 
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Theorem 1.6. The identity (21) is the “shadow” of a similar summation 


identity for the Eisenstein series. That is, for n > 1 we have 


GQ n +2^z') — 


n 


1 (4n + l)! 

6n+ 1 ' (2n)! 2 ^ 


/ 2 n \ 

T—— rG2n+2k{z)G^ n -2k+2{z). (23) 

V2n+2fc-l/ 


Let us now explain the connection between our results on co(s) and the 
identity (20). Remarkably, we will show that the relation © identifying 
the “trivial zeros” s = —1, —2,... of uj(s) is equivalent to that identity. It 
was precisely our study of the values w(—n) which prompted the discovery of 


(20) as a numerical observation for small values of n. We were then able to 


prove the identity (and the stronger Theorem 1.6) using separate algebraic 


methods and deduce (17), and also learned about the earlier work of Agoh 
and Dilcher. The same reasoning could work in reverse to deduce (20) in 


yet another way if a more direct proof of © could be found (e.g., based 
on a functional equation satisfied by cj(s)), but we do not currently know of 
such a proof. 

To conclude this discussion, let us comment on the connection to the the¬ 
ory of modular forms. For k > 0, let denote the vector space of modular 
forms of weight 2k for the modular group SL( 2, Z). The space M 2 k is known 
to be of dimension [k/6 J if k = 1 (mod 6), or \k/ 6j +1 otherwise. The 
Eisenstein series Gq u+ 2 , and the products G 2 n + 2 k(z)G 4n - 2 k+ 2 {z) appearing 


on the right-hand side of (23), are all elements of M$ n+ 2 , and, curiously, 


the number of distinct terms on the right-hand side of (23) (after group¬ 


ing together the equal pairs G 2 n +2k{z)G 4 n -2k+2{z), G 4 n -2k+2{z)G 2 n +2k{z) 
obtained from symmetric indices k and n + 1 — k) is precisely equal to 
dimM 6 n+ 2 ; check this separately for both even and odd values of n. Thus, 
the theory of modular forms seems like the proper context in which to con¬ 


sider (23) and similar summation identities. It is natural to wonder whether 


there is a deeper connection between (23) and the function cj(s); whether 
other connections of this type between modular forms and zeta functions can 


be found; and whether summation identities such as (23) can be classified 


and better understood. See also the list of open problems in Section 12 


1.4 Methods, organization, and the companion Mathematica 
package 

The remainder of the paper is organized in roughly three parts: the first part 
consists of Section [2] and [3j in which we give a new proof of the summation 
identity (20), and then show how a modification of our argument yields a 
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proof of Theorem 1.6 Our method is an extension of an algebraic technique 
introduced by Zagier [50] . that enables us to reduce the identities to a family 
of more conventional binomial summation identities. These identities are 
proved using Zeilberger’s algorithm and the method of Wilf-Zeilberger pairs. 

In the second part of the paper, comprising Sections [4] and [5j we study 
the properties of the function w(s), and prove Theorems 00 


and 1.4 


The analysis relies heavily on an integral representation for cu(s) derived by 
Matsumoto [25] . 

In the final part, which consists of Sections mm we apply the results 
concerning cj(s) to the study of the asymptotic behavior of the sequence r(n). 
We combine Mellin transform methods with a saddle point analysis (framed 
in probabilistic language as a local central limit theorem), culminating in 
the proof of Theorem mi A key number-theoretic difficulty we encounter 
in the course of the proof is that of proving effective bounds for the decay 
of a certain generating function away from the saddle point; we derive the 
necessary bounds by using the modular transformation properties of the 
Jacobi theta functions. We conclude with some final comments and open 
problems in Section |T2j 

A few of the proofs in the paper involve tedious computations of a purely 
algebraic nature, which can be checked by hand without much difficulty but 
are most easily verified using a computer algebra package. We prepared a 
companion Mathematica package that is available to download online 
and demonstrates the correctness of these computations. 
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2 Bernoulli numbers 


In this section we give a new proof of the summation identity (20) using a 
completely different method than that used by Agoh and Dilcher [2]. In the 
next section we show how a modification of our argument actually proves the 


stronger identity (23) involving the Eisenstein series, which contains (20) as 


a specialization at z = zoo. 
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Note that we gave three equivalent formulations (20)—(22) of the summa¬ 


tion identity. While (20) seems like the most natural number-theoretic for¬ 


mulation of the identity (and was the version proved by Agoh and Dilcher), 
the other two formulations will also play important roles, in that our proof 


technique will actually imply the version (21), whereas we will later make use 


of its “dual” version (22) (the two being related via the functional equation 
for the Riemann zeta function) in our analysis of the function uj(s). 

Our proof is based on adapting a technique introduced by Zagier [50] as 
a way of proving the more elementary recurrence relation 


71—1 


C( 2 ^) = 2n j r \ Y C( 2fc )C(2n - 2k) (n> 2), 


k =i 


a well-known relation due to Euler. (If one assumes the connection between 
the even integer zeta values to the Bernoulli numbers, this is the recurrence 
B 2 n = 2 n+i J2k=i ilk) B 2kB2n-2k, which has a simple proof using gener¬ 
ating functions; see [8], Proposition 1.15].) Before we proceed with the proof 
in the most general case, let us first illustrate the use of Zagier’s technique 


to prove the case n = 1 of (21), which in that case becomes 


C(8) = fC(4) 5 


Define the function f{p,q) = ^5 + 7 ^ 3 - + One may easily verify 

that f(p , q ) satisfies the algebraic identity 

6 

f(p,q) - f(p + q,q) - f{p,p + q ) = yrr- 

7jrq^ 

Now sum both sides of this equation over all 6 N, The summation of 
the right-hand side gives 6£(4) 2 /7, and for the left-hand side we get 


5^/(p,?)- Y f(p + ^,q)~ Y /(p»p + «) 

p,q=l p,q= 1 p,q =1 


Y - Y - Y /(*■>*) = £/(*■>' 


r,s=l r>s>l s>r>l 

2 3 2\ ^ 1 

7 r =1 


r=l 


For general to, we now try to generalize the idea by looking for a pair of 
functions f n (p, q), g n {p, q) = fn{p , q) ~ fn{p + q,q)~ fn{p,p+ q), such that 
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f n (p , q ) is of the form 


2n+l 

fn{p,q) = E p 2n+j q An-j+2 ’ ( 24 ) 

J=1 

for some as-yet-undetermined coefficients (a' n j) 2 ”'| 1 , and such that g n (p,q) 
miraculously simplifies to be of the form 

9n(Pi q) — E 2n+2k a An-2k+2 ’ ( 2 ^ 

fc=l 1 1 

for some coefficients (/3 n fc)fc=i- Assuming that this somewhat daring ansatz 
turns out to be successful, summing g n (p- q) over all p, q G N then yields on 
the one hand 


n 

£ 

fc=i 


| E p - 2 "- 2 fc g - 4 - +2fc - 2 

V.P.9= 1 


E/?n,fcC(2n+2A:)C(4n-2A:+2), (26) 

fe=i 


and on the other hand 


E fn(p,q)~ E /^(p + ^9)- E /-(P’P + '?) 

p,g=i p,9=i p,g=i 


E - E - E /-(as) = E/-( r >' 

r,s=l r>s>l r>s>l J r =1 

( 2n+l \ 

E J' j C(6n + 2), 


(27) 


whereupon we obtain, after also adding a normalization condition )Th a n j = 
1 , the identity 


n 

C(6 n + 2) = E ^n,fcC(2n + 2fc)C(4n - 2fc + 2). (28) 

fc=i 

The difference with Zagier’s original application of the technique is that 
in his case the undetermined coefficients were trivial (see the bottom of 
p. 498 of his paper [50]), whereas here they are not. As the next lemma 
shows, finding them reduces to the solution of a somewhat tricky system of 
linear equations. 
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Lemma 2.1. If f n (p, q) is defined by (24), then the function g n (p,q) = 
fn(,P, q) ~ fn(p + Qi q) ~ fn(p , P + q) can 6e expressed in the form p5] ) i/ and 
only if the coefficients (anj) 2 ^ 1 , (fin,k)k=i satisfy the equations 


2n+l 

E 

i=i 


a 


n,J 


6 n 

m — 4n + j 


2 n + j — 2 
m. — 4n + j 




fc=i 


4n — j 
rn — 6 n 


6 n 

m — 4n + 2fc 


/or m = 2n, 2n + 1,..., lOn — 2. 


(29) 


Proof. In (25), substitute the definition of g n (p,q ) and multiply both sides 
by p 6n g 6n (p + g) 6 ". The identity becomes 


2n+l 

E 

3= 1 


^ a n j p An -jq 2n+ i- 2 (p + q) 6n - p An ^q 6n (p + q) 2n+ *~ 2 


- p 6n q 2n+j ~ 2 (p + qf n -l = Y, fin,kP An - 2 k q 2n+2k ~ 2 (p + qf n , 

) k=1 

an equation in which both sides are homogeneous polynomials in p, q of 
degree 12n —2. Extracting a common factor of q 12n ~ 2 and denoting x = p/q, 
the identity is clearly seen to be equivalent to the single-variable polynomial 
identity 


2n+l 

E 

3= 1 


a 


n,3 


x 4n ~ j (l + x) 6n - x 4n ~ j (l + x) 2n+j ~ 2 - x 6n (l + x) 4n ~ j 


E&. 

k= 1 


kX 


.An—'2k 


(l + x ) en . 


Expanding both sides of this relation in powers of x and equating coefficients 
of each monomial x m gives precisely the equations (29). The power m on 
the right-hand side ranges from 2 n to lOn — 2; on the left-hand side it ranges 
from 2n — 1 to lOn — 1, but it is easy to verify that the coefficients of x 2n ~ l 
and x 10 "- 1 vanish. □ 


Having derived the system of equations (29), we used a computer to 
solve them for small values of n. This led us to guess based on numerical 
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evidence that 


a n,j ~ 


fin,k 


(4n+l)! (A) 


6n+l (2n!) 2 

2 n 


2 (4n + l)! ( 2 ^- 1 ) 

6n + r (2 n!) 2 ( 2n+ 6 2 Vr) 


= 2a 


n,2k' 


(30) 

(31) 


In our situation we happened to already know what the coefficients 


should be, since we knew in advance the form of the identity (28) to be 


proved—see Section [3] for the reasoning that led to its discovery. However, 
it is worth noting that the equations can be solved even with undetermined 
values for these coefficients, which can be useful for discovering additional 
identities; see Section [12] where we give two additional examples of identities 
discovered using this approach. 


Having guessed the formulas (30)—(31), our final task will be to verify 


that the guesses are correct. We start by checking the normalization condi¬ 
tion. 


Lemma 2.2. If the numbers a n j are defined by (30), then for all n > 1 we 
have that Yl“i= l" 1 a n,j = 1 - 


Proof. Start by observing that the definition (30) of the a n ,• ’s as an inter¬ 


esting triangle of integers also makes sense when n is a half-integer. It turns 
out that the identity Ylj a n,j = 1 is still correct in that case, and this is 
both easier and more natural to prove. That is, we now wish to show that 


S(n):='£F(n,k) = 1 (32) 

k =0 


for all integer n > 0, where we define 


1^(71, k ) Q ! n/2,fc+l 


1 (2n + 1)! $ 

3n + l (n!) 2 Q k )' 


This can be proved using the method of Wilf-Zeilberger pairs [57j. Define 
the “proof certificate” function 


R(n, k ) = 


k(2n -k + l)(lln 2 + 2 k 2 - 5 nk + 27 n -6 k + 16) 
3(n + l)(3n + 2)(3n + 4)(n - k + 1) 


which was found with the help of the Mathematica package fastZeil [37], 
|38j, a software implementation of Zeilberger’s algorithm [53] (see also 
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Let G(n,k ) = F(n,k)R(n,k). Then one may verify by direct computation 
that 

F{n + 1, k) — F{n , k ) = G(n , k + 1) — G(n, k) (33) 

(divide both sides by F(n , fc) to turn the relation into an equation between 
two rational functions, or refer to the companion package [33j for an auto¬ 


mated verification). Summing both sides of (33) over all integer k shows 
that S(n+ 1) = S(n). Since 5(0) = 1, (32) follows by induction. □ 


Lemma 2.3. The coefficients a n j and j3 n) k as defined by (30) —(31) satisfy 
the system of equations (29). 

Proof. Since j3 n) k = 2a n ^k, it makes sense to rearrange the terms between 
the two sides of the equation (29), arriving at the equivalent relation 

4n — j 


2n+l 

£ 

3 = 1 


(-D 


i-i 


6 n 

m — An + j 


a n,j ~ 


2n+l 

£ 

3 = 1 


2 n + j — 2 
m — An + j 


+ 


m — 6n 


a. 


n,J- 


( 34 ) 

It is now worth noting for aesthetic reasons, although not strictly necessary 
for our purposes, that, as occurred in the proof of Lemma 2.2 above, a 
version of this identity also holds when n is a half-integer, although to get the 
correct identity it is necessary to insert in two places a small sign correction 
of ( — l) 2n . Making this correction, and substituting a U) k into the equation, 
with the irrelevant scaling factor (which does not depend on j ) 

removed from both sides, the identity we need to prove becomes 


2n+l 

£<-u 

3 = 1 


2 n+j 

(A") 

/ 6 n ^ 

\2n-\-j— l) 

2n+l 

(2n\ 

V 

\j-l) 


( 6 n 

3 = 1 

\2n-\-j— 1 


6 n 

m — An + j 

2n + j — 2 
m — An + j 


+ (-l) 


2 n 


An — j 
m — 6n 


(35) 


Now replace n by n/2, and perform two additional simplifications of a purely 
cosmetic nature, namely shifting the summation index j by 1 and replacing 
the parameter mbym + n-1, This brings (35) to the form 

Of 3n 

\m + j — n 

(?) 


£(-d” + L 3 „, 

\n+7/ 


3=0 


n+jJ 

n 

£ 


( 3n .) 

=0 \n+j) 


n + j-1 
m + j — n 


+ (-l) ? 


2 n — j — 1 
in — 2n — 1 


(36) 
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where n is assumed to be a positive integer and m ranges over the values 
0,1, 2,..., 4n. 

Our next step is to exploit a subtle symmetry of the family of identities 


(36) in order to obtain one final simplification. The key observation is that 


the binomial coefficient ( m _ 4 7 l n +j) on the left-hand side is invariant under 
the substitution j —> n — j, m -+ 4n — m. On the right-hand side, the same 
substitution transforms the two binomial coefficients and (^i 2n ^\) 

into each other. Note further that when 0 < m < 2n — 1, only the first of 
those two coefficients can be nonzero (under the constraint that 0 < j < n), 
and when 2n + 1 < m < 4n, only the second one can be nonzero. When 
m = 2n both of them are zero, and in this case it is easy to verify that the 
left-hand side is a sum with alternating signs of the binomial coefficients (”) 
and therefore also zero. Thus, by symmetry it is enough to verify ( |36[ ) for 
0 < m < 2n — 1, and for that range of values it becomes the relation 


E 

3=0 


( n ) / 3 n \ 

(_1 ^n+j}3' \m+j—n) 


( 3r M 

\n-\-j) 


E 

3=0 


/n\ / n+j—1\ 
\j J \m+j—n) 

{ 3n .) 


(37) 


Following this preparation, the well-known algorithmic methods for proving 
hypergeometric summation identities—specifically, Zeilberger’s algorithm— 
once again come to our aid. Denote 

/n\ / 3n \ n 

Fi(n,m,j) = (-l) n+J J , Si{n,m) = ^F 1 (n,m,j), 

\n+j) 


F 2 (n,m,j) = 


(n\ /n+j—1\ 
\j ) \m+j—n) 

pM ’ 


3=0 

n 


S 2 (n,m) = y^F 2 (n,m,j), 

3=0 


and define the additional auxiliary functions 


U(n, m ) = (m + 2)(2n — m — 1), 

V ( n , m) = —2 nr? — 5n 2 + 8 mn + 9 n — 4m — 2, 
W(n, m ) = (m — 4n)(2n — m — 1), 
j(3n + l)(2n- j + 1) 


n — m — j — 1 
j(m — 2 ra + l)(2n — j + 1) 


Ri(n,m,j) = 

R 2 (n,m,j)= . 

n — m — j — 1 

Gi(n,m,j) = R 1 (n,m,j)F 1 (n,m,j), 
G 2 (n, m, j) = R 2 (n,m,j)F 2 (n,m,j). 
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The functions U (n, m), V (n, m), fT(n, m), Ri(n, m, j), and R. 2 (n, m, j) were 
found using Zeilberger’s algorithm, and one may verify by a routine calcula¬ 
tion (see the companion package [T3j) that they satisfy the pair of algebraic 
identities 


U (n, m)F a (n, m + 2, j) + V ( n , m)F a (n, m + 1, j) + W(n, m)F a (n, m,j) 

= G a (n,m,j + l)-G a (n,m,j) (a = 1 , 2 ), (38) 

where, importantly, the polynomials U(n,m),V(n,m),W(n,m) do not de¬ 
pend on j or a. Summing both sides of (38) over all integer j results in 
a telescoping sum on the right-hand side, which (since G a (n,m,j) = 0 if 
j < 0 or j > n) is equal to 0. Thus, we get the relations 

U(n , m)S a (n, m + 2) + V(n, m)S a (n, m + 1) + W(n, m)S a {n , m) = 0 (39) 

for a = 1 , 2 . That is, we have shown that the sums Si(n, m), 62 ( 71 , m) 
satisfy the same second-order linear recurrence relation with respect to the 
variable m. 

Finally, note that U (n, m) ^ 0 for 0 < m < 2n — 2, so for those values of 
m we can solve the recurrence ( |39[ ) for S a (n, m + 2 ), expressing it in terms of 
S a (n, m + 1) and S a (n, m). We therefore see that it is enough to verify that 
S i(n, 0 ) = 52(71,0) and <Si(n, 1 ) = ^2(77, 1), after which the general equality 
Si(n,m) = 52 ( 77 , 777 ) (which is precisely the relation (37) we are trying to 
prove) follows for all 0 < 777 < 2 t 7 by induction on m . The verification of the 
base cases 777 = 0,1 is trivial and is left as an exercise. □ 


Combining Lemmas 2.1, 2.2 and 2.3 establishes the summation identity 
(28) with the coefficients f3 nt k as defined in (31), and therefore finishes the 
proof of (p0 ). 


3 Eisenstein series 


The key idea in our proof of the Bernoulli number summation identity (20) 


in the last section was to find a pair of functions f n (p , q), g n (p, q) = f n (p, q) ~ 
fnip+q, q)~fn(p,P+q) which are both homogeneous polynomials in p~ 1 ,q ~ 1 
of degree 6?7 + 2, and such that f n (p , q) contains only a specific set of even 
powers of p -1 and q -1 . Zagier, who introduced this method in [50j to prove 
a more standard identity, already made the more general observation (see 
section 8 of his paper) that any summation identity proved using the same 
technique automatically “lifts” to a summation identity for the Eisenstein 
series, by a modification of the argument. To see how this works in our 
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setting, instead of summing g n (p- q) over all pairs of positive integers p, q as 
we did in (26), (27), now sum instead over all pairs of complex numbers p,q 
ranging in the “half-lattice” 

A_|_(z) = {p = mz + n : rn > 0 or [rn = 0 and n > 0]} 

(where z is a fixed complex number in the upper half-plane H, which will 
shortly become the argument of the Eisenstein series). This summation 
yields, on the one hand (analogously to (26)), the expression 

G , 2n+2fc(") G < 4 n _2fc+2(^) 


E« 

fc=l 


’n,k 



p -2n-2k q -4 n+ 2k-2 | = ^ ^ 

k =1 


On the other hand, for two elements r = mz + n and s = m!z + n' of the 
half-lattice A_|.(z), denote the order relation r s to mean that m > m' 
or [m = rn' and n > n '}, and note that r y s if and only if r = s + q for 
some q € A+(z). Then the summation above can also be written, in analogy 
with (27), as 

OO OO OO 

^2 fn(p,q)~ fn(p + q,q)~ ^2 f^{p,p + q) 

p,q£A + (z) p,q£A+(z) p,q£A+(z) 


\fn(r,s)= ^2 fn(r,r) 

r£A + (z) 


- 


OO 

E - 

E 

- E 

r,s£A + (z) 

r,seA + (z) 

r,sEA+(2:) 


rys 

syr 

2n+l \ 

22 an ’jJ 

GrQ n - j-2(^) 

G 611 + 2 ( 2 ) 

2 

2 


Thus we obtain precisely the identity (23), proving Theorem 1.6 


□ 


In connection with the proof above, we remark that Zagier’s observation 
about the lifting of summation relations to Eisenstein series was made in 
the context of considering general summation relations of the form 


G2k — 22 C jG2jG2k-2j 

2 <j<k 


(40) 


expressing a given Eisenstein series in terms of the Eisenstein series with 
lower index. It is apparent from the formula for dim M 2 k mentioned in 
Subsection 1.3 that for any fixed k there must be a certain number of linearly 
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independent relations of this kind (the number of such relations increases 
linearly with k), and one can use standard algorithms from linear algebra 
and the theory of modular forms to find all such relations. From this point 
of view, the existence of summation relations of the form (40)—again, for 
specific values of k —is not especially surprising. 

What is not necessarily apparent, and it is not clear to us from the dis¬ 
cussion in m whether it was envisioned by Zagier, was that his technique 
can also be used to prove interesting families of summation relations (other 
than the single relatively trivial example discussed in section 1 of his pa¬ 
per) that hold for many values of k and have the kind of elegant structure 
that our identity (23) possesses. The discovery that this is in fact the case 
makes Zagier’s technique appear a lot more exciting, and opens the way 
to discovering and proving—possibly in an algorithmic fashion—additional 
families of identities; see Section [12] for a further brief discussion of these 
possibilities. 


4 The function tu(s) and its analytic continuation 


The goal of this section is to prove Theorems |1.2| and |1.3[ Our approach is 
based to a large extent on an integral representation derived by Matsumoto 


for the more general Dirichlet series Cmt(si, S 2 , S 3 ) defined in (14) (see equa¬ 
tion (5.3) in [25]). 


The region of convergence: proof of Theorem |1.2| (1). We start by 
giving a simple argument to prove the claim of Theorem |1.2| (1) about the 
precise region of convergence of the Dirichlet series defining w(s). Observe 
that 


^ U k U + k )) s 


■( 


E 


1 


\j>k> 1 


(jk(j + k)Y 


1 


+ E 


1 


2 V — , v 

j>k>i j 2s k s (1 + j) 


k>j> 1 

+ 2" s C(3s). 


(jk(j + k)f 


+ E 

k=j>l. 


1 


(jk(j + k)) s 


Here, the series for 2 S £(3s) converges if and only if Re(s) >1/3, and the 
summation over j > k > 1 converges absolutely if and only if the simpler 
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series 


1 


(41) 


oo 1 

y _L 

I j^ s k s 

j>k >1 


E 

_7>fc>l 


j 2a k° 


(where a = Re(s)) converges absolutely, since on the range of summation 
the multiplicative factor (l + y) s is bounded between 2~ s and 1. It is now 
trivial to check that the series in (41) converges absolutely and uniformly 
on compacts (which ensures that the sum is a holomorphic function) if and 
only if Re(s) >2/3, for example by comparison with the integral 

1 , , 1 / 


n 


x 2a y c 


dy dx = 


„1 —3er 


(J — 1 


dx. 


□ 


Analytic continuation: proof of Theorem |1.2| (2). Next, we proceed 
with the analytic continuation of uj(s) to a meromorphic function on C. We 
recall Matsumoto’s ingenious method to achieve this (which he applied in a 
more general context), but add a few more details and motivation that are 
missing from his analysis and that we feel add some important insight in 
view of our newer discoveries about the significance of a;(s). Start with a 
standard contour integral formula from the class of Mellin-Barnes integrals 
[36, eq. (3.3.9), p. 91], namely 

h(s)(l + A) -s = ^ / h(a + z)T(-z)X z dz , (42) 

(a) 

where (a) denotes the vertical contour from a—ioo to a+ioo (with i = \/—l, 
here and throughout), and where we assume that A G C \ (—oo,0] and 
—Re(s) < a < 0. Setting A = k/j, multiplying by j~ 2s k~ s and summing 
over j, k G N gives that 



= 2 J F(s +z)T(-z)((2s + z)((s - z)dz, (43) 

(a) 
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provided the infinite summations converge and that the interchange of the 
summation and integration operations can be justified. Note that the two 
zeta series within the integrand can be summed precisely when Re(2s + 2 ;) > 
1 and Re(s — z) > 1, that is, when 1 — 2a < a < a — 1 (where a = Re(s) 
as before). Together with the constraints —a < a < 0, we see that a and a 
must satisfy 

max(— a, 1 — 2a) < a < min(0, a — 1). (44) 

As illustrated in FigureQa), a feasible value of a satisfying these constraints 
exists if and only if cx > 2/3, which is consistent with our earlier observations 
on the region of convergence of the Dirichlet series defining cu(s). If a > 2/3, 
and a valid value of a is chosen (for example, the value a = —1/3 works), the 
interchange of the summation and integration operations is easy to justify 
using standard facts about the exponential decay of the gamma function, 
and the (at most) polynomial growth of the zeta function, along vertical 
lines, namely that for z = x + iy, \y\ > 1, we have that 


\r(x + iy)\ = C>(\y\ A e-^A, 

(45) 

\({x + iy)\ = 0(\y\ A ). 

(46) 


where A is a constant that depends only on x, and where both A and the 
constant implicit in the big-O notation are uniform as x ranges in compact 
subsets of M; see [36/ Section 2.1.3] and [M, Section 6.4] for the proofs of 
these relations. 


Summarizing the above discussion, we conclude that the integral in (43) 
is indeed a valid representation for T(s)a;(s). We remark also that the 


bounds (45) and (46) together with standard facts from complex analy¬ 


sis will also provide easy justifications for all subsequent manipulations of 
contour integrals in this paper (e.g., shifting the contour and differentiation 
under the integral). We assume the reader is familiar with such techniques 
and will mention these justifications only in passing without going into de¬ 
tail. 

The next step is to analytically continue the integral to a larger region. 
This is done by shifting the contour of integration and taking into account 
the residues. Note that the integrand has poles (as a function of z with s 
fixed) at: 

• z = s — 1 because of the factor ((s — z); 

• z = 1 — 2s because of the factor £(2s + z); 


21 





• z = 0,1,2,... because of the factor T(— z)\ 

• z = — s, — s — 1, — s — 2,... because of the factor T(s + z). 

This means that we should look for a value of a for the new integration 
contour (a) for which a = Re(s) is allowed to range over a wide range of 
positive and negative values without intersecting the lines in the a-a plane 
defined by the equations 


a = a — 1 , 


(47) 

a = 1 — 2ct, 


(48) 

a = n 

(n = 1,2,...), 

(49) 

a = —o — n 

(n = 0, !,...)• 

(50) 


The situation can be visualized graphically in terms of the “singularity di¬ 
agram” shown in Figure Qb), which shows the pairs (a, a) that cause the 
integral to pass through a singularity point. 


We are ready to shift the integration contour in the representation (43). 
Simplifying Matsumoto’s argument slightly, we assume for convenience that 
initially | < a < 1, so that we can avoid having to consider the double 
poles that result when two of the singularity lines intersect. The diagram 
in Figure Qb) shows that it is a good idea to shift the contour (a) to the 
right (moving up in the coordinate system of the diagram) to (a!) where 
a! = M — 1/2 and M is some large integer. As the diagram shows, this 
involves skipping over the pole at z = s — 1, which is a pole with residue 
—r(2s — l)r(l — s)C(3s — 1); and the pole at k for each k = 0, 1, 2, ..., M — 1, 
which has residue —T(s + &)£(2s + k)((s — k ). Thus, by the residue 


k\ 


theorem we have (still assuming g < a < 1) that 


r(s)w(s) = r(2s - i)r(i - s)c(3s -1) 

/ i\fc 

+ ^2 ^ r(s + /c)C(2 s + k)((s — k) 

k =0 

+ [ r(s + z)r(—z)£(2s + z)£(s - z) dz. (51) 

2m J 

(M—1/2) 


Furthermore, Figure[l|b) shows graphically, and the equations (47)—(50) can 
easily be used to show algebraically, that the integrand in (51) encounters 
no singularities as s ranges in the strip 


3 

4 


M 

Y 


< a < M + 


(52) 
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Figure 1: (a) The shaded area in the diagram represents the feasible region 
of pairs (a, a) for which the constraints (44) are satisfied and therefore the 
representation (43) is valid, (b) The lines in this diagram represent pairs of 
(a, a) where the integrand in (43) will run into a singularity. To analytically 
continue cu(s), we start at (cr,—1/3) where a is initially allowed to vary in 
the interval (|, 1), (the triangular shaded region) and move up along the a 
axis (which corresponds to shifting the contour of integration to the right) 
to a region of the diagram (the trapezoidal shaded area) where a can be 
varied across a much wider interval without running into singularities. 


It follows (using 


-(46)) that the integral defines a function which is 
Since M was arbitrarily large, we have 


holomorphic in the same strip 
therefore succeeded in extending ca(s) to a meromorphic function on the 
entire complex plane, proving Theorem |1.2| (2). □ 

Analysis of the poles and values at nonnegative integers: proof of 
Theorem |1.2| (3), Theorem |1.3[ We now extend Matsumoto’s analysis 
further by examining the singularities of w(s) and its values at nonnegative 
integers. First, consider r(s)w(s), whose analytic continuation is given by 


the right-hand side of (51). Here is a list of the poles of r(s)cu(s) contributed 


from each multiplicative factor in each of the additive terms, along with their 
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residues. 


1. The poles of the term T(2s — 1)T(1 — s)£(3s — 1) are: 

• A simple pole at s = § because of the factor £(3s — 1). The 
residue is |T Q) 2 . 

• Potential simple poles at s = 0, — — 1, — — 2,... because of 

the factor T(2s — 1). The residue at s = for each integer 

k > 0 is ^ 2 -fci r (^ 2 ^) C ( 1 ~ 2 3fc ) • Note that for k = 3, 7,11,15,... 
the residue is 0 because of the trivial zeros of ((s) at the negative 
even integers, so the poles of the term T(2s — 1)T(1 — s)£(3s — 1) 
at these values are in fact removable singularities. 

• Simple poles at s = 1, 2, 3,... because of the factor T(1 — s). The 
residue at s = k + 1 is equal to -—— T(2k + 1)£(3 k + 2) = 
< '~ 1 ^ k\ < ' 2fc ^ C(3fc + 2) for each integer k > 0. 

2. The poles of T(s + k)((2s + k)((s — k) for k = 0,1,..., M — 1 are: 

• A simple pole at s = k + 1 because of the factor ((s — k). The 
residue is T(2 k + 1)£(3 k + 2) = ^ 1 \.P fc ^ C(3fc + 2). 

• A simple pole at s = because of the factor £(2s + k). The 
residue is (=^r(^)C(^). 

• Potential simple poles at s = —k, —k — 1, — k — 2,... because of 
the factor r(s + k). The residue at s = —k — j for each integer 
j > 0 is ~ (-j? J ((—k — 2j)C(—2k—j). As above, some of them 
end up being removable singularities because of cancellation with 
a trivial zero of C( s )- 

Adding up the contributions from the different terms (assuming that 
M is arbitrarily large so we can ignore the edge effects), we see that the 
poles of r(s)w(s) at s = 1,2,3,... have residue 0 and are in fact removable 
singularities (which could have been predicted in advance, as we already 
knew that T(s)cu(s) is analytic at those points). The remaining poles consist 
of a “special” pole at s = | with residue (^) , and a series of (potential) 
poles at half-integer points of the form s = for k = 0,1, 2,..., where for 
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an even value k = 2m the residue is 

(_l)2m /I _ 6 m 

(2m)! v 2>s 

1 • 3 • 5 • ... • (2m A - 6m 

(.2 ) C 


(2?n)!2 m 

and for k = 2m + 1 an odd integer, the residue is 

f _j\2m+l 

Pm := 75-^TTT r ( m + l)C(-3m - 1) 


1 

4 m rn! 


v^C 


1 — 6m 
2 


(2m + 1 )! 


+ 


f-1) 


m\ 


3=0 


( j )C(-ra- j)C(-2m + j) 


m ' f C(-3m - 1) + ^ j ^ ) C(-"i - j)C(-2m + j). 


(2m + 1)! 


m! 


j=o 


Note that po = C(0) 2 — C( — 1) = (— |)~ + A = A However, as we discovered 
to our surprise, the remaining values p\,p 2 ,--- are all zero! For odd values 
of m this is trivially true since ((—3m — 1) is a trivial zero of the Riemann 
zeta function, and in the summation over j. for each j one of the factors 
(depending on the parity of j ) is also a trivial zeta zero. For even values 
of m, denoting m = 2 p, we see that the claim reduces to the summation 
identity 


((—6p — 1) = (4p + l) 


4 p 

2 p 


2 v 

E 

j =o 


2 p 


((-2p-j)((-4p + j), 


which is precisely the version (22) of our Bernoulli number summation iden¬ 


tity (note that the summation over j can be taken only over odd values, 


again because of the trivial zeta zeros). Thus, (20) and its reformulations 


(21), (22) turn out to be equivalent to the claim that r(s)cj(s) has no poles at 
the negative integers. Since we, and Agoh and Dilcher before us, proved the 
former result using other methods, the latter claim about r(s)w(s) follows. 

We summarize our findings about r(s)w(s) in the following theorem, 
which will play an important role in our asymptotic analysis of the repre¬ 
sentation enumeration function r(n). 


Theorem 4.1. The meromorphic function T(s)w(s) possesses a simple pole 
at s = | with residue |r Q) ; a simple pole at s = 0 with residue | ; and for 
each k = 0,1, 2,... a simple pole at s = \ — k with residue ^( 1 ^ 6fc ) ■ 
It has no other singularity points. 
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From Theorem 


4.1 we can now easily infer part (3) of Theorem 1.2, and 


both claims (16), (17) of Theorem 1.3 This is based on two observations: 


first, for a complex number so which is not a pole of the gamma function, 
we have 

Res s= <j 0 (w(s)) = y Res s=So (r(s)w(s)). 

Second, at the poles so = —n- of T(s), n > 0, we have 

(r(a)w(a)) 


w(so) = 


Res s=so 


Res s=So 


= (—l) n n! Res s= _ n (r(s)o;(s)) 


(where the equation holds in the usual complex-analytic sense, that is, after 
removing the removable singularities). Applying the first observation to 
so = — 5,..., and the second observation to so = 0,-1,—2,..., 

gives the claims; for example, for sq = | we get that 


Res s=2/3 (M>o) = ir(i)Vr(l), 


which is equal to 9 J^ T Q ) 3 by the reflection formula r(x)T(l— x) = 

We leave to the reader to fill in the details for the remaining poles of 

r (s)w(s). 


We have therefore proved parts (l)-(3) of Theorem 1.2, proved Theo¬ 


rem 


1.3 and established the additional interesting fact that (16) implies, 


and is implied by, the summation identity (20). Furthermore, note that the 


evaluation w(0) = 3 is directly related to the fact that the identity (20) does 


not hold for n = 0. 


□ 


Polynomial bound along vertical strips: proof of Theorem |1.2| (4). 

To conclude our analysis of w(s), it remains to prove part (4) of Theo- 
namely the claim that |ca(cr + it) \ grows at most polynomially in t 
00 with a ranging in a bounded interval I = [a, b]. This fact was 


1.2 


rem 
as |t| 

mentioned briefly in |26| (in the more general context of the multivariate 
Mordell-Tornheim zeta function) and is a fairly straightforward consequence 
of the representation ©: dividing that representation (with M taken large 
enough so that the representation is valid for a £ I) by T(s) and using the 
standard relation r(u; + 1) = wT(w), we get that 
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U)ts = 


r(2s- i)r(i-s)C(3s- i) 


M—l 

+ £ 

k=0 

1 


r(s) 

(-i)* 

k\ 


s(s + 1)... (s + k — l)£(2s + k)C(s — k ) 


+ 


2iri 


r(s + z)r(-z)c(2s + -z)£(s - z) 
r(s) 


dz. 


( 53 ) 


(M—1/2) 

Here, all the summands except possibly the integral can be immediately 


seen using {45)—(46) to satisfy such a polynomial bound. For the integral, 


writing s = a + it and z = x + iy, we use (45)—(46) again to see that it can 
be bounded as 


(M—1/2) 


T(s + z)T(-z)((2s + z)C(s - z) 

fw 


dz 


7T j . 

< e* 1 


/ OO 

e -2 It+s/lg-f hi x |p 0 ly(t, y)| dy , (54) 

-OO 


where Poly(t, y) denotes a polynomial factor (with powers and coefficients 
that are uniform as a ranges over compact intervals). However, a quick side 
calculation, which we omit, shows that for any integer n > 0 the function 


«i->e 


r 


|mng-M e -h+«l du 


is a polynomial in u. Applying this to the right-hand side of (54) gives the 
desired polynomial bound, and finishes the proof of Theorem 1.2 □ 


5 The formula for u/(0): proof of Theorem 1.4 


To prove the formula (18) for o/(0), start with the integral representation 


531). In the case M = 2, according to (52) the representation is valid when 


— -j < Re(s) < §. In particular, we have in a neighborhood of s = 0 that 


uj(s) = 


T(2s — 1)T(1 — s)£(3s — 1) 


T(s) 


+ C( s )C(2s) — «C( S — l)C(2s + 1) 


+ 


2tti 


(3/2) 


T(-z)T(.s + ^)C(2s + z)((s - z) 

W) 


dz, 


(55) 
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where in the case of the removable singularity at s = 0 the equality is under¬ 
stood in an appropriate limiting sense. Now inspect the Taylor expansion 
of each of the summands up to first order, using standard facts about the 
Taylor expansions of the gamma and zeta functions at s = —1, 0,1. For the 
first summand, we have 

T(2s-l)T(l-a)C(3a-l) 

_ (1 + 7 s + 0( g2 )) (-^ + (7 - 1) + Q( 5 )) (~T 2 + 3C ; (-i)s + Q( fi2 )) 

7 -l + 0{s) 

= ^ + (^ _ ^ c ' (_1) ) s + 0(s2) - 


Similarly, for the second summand, 

C(s)C(2s) = ^ log(27r)s + 0(s 2 )^j - log( 2 vr)s + 0(s 2 )^J 

= \+ ^l°g(2vr)^ s + 0(s 2 ). 
and for the third summand, 

-<(s - l)C(2a + 1) = -a ^ + 7 + O(s)^ + C'(—1)« + 0 (s 2 )^ 

= ^ + (^-^' ( - 1 >) >+0<A 
The coefficient of s in the sum of the above three expansions is 

12 + T2 7 + I log ( 27r ) - 2 ^( -1 )’ ( 56 ) 


which accounts for all the terms on the right-hand side of (18) except for 
the integral. 

It remains to consider the integral in (55). Define 


A(s,z) = 


T(—z)T(s + z)((2s + z) C(s - z) 
_ . 


Differentiating A(s,z) with respect to s and using the Taylor expansions 


F(a) 

F'(s) 

T(s) 


- — 7 + 0(s), 
s 

-- +7 + 0(s), 

s 
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a short computation, which we omit, shows that 


d_ 

ds s=o 


A( s , 2 ) = r(,)rH)c(z)(H 


which, using the functional equations T(u; +1) = and r(’:)r(l — z) = 

7 t/ sin(7r^), simplifies further to 


7T-I A(S,Z) =- — ^((z)((~z). 

as |s=o zsm{TTz) 


Thus, differentiating the integral in (55) and setting s = 0 gives 


d_ 

ds | s =o 

1 

~2i 


1 

2m 


T(—z)T(s + z)C{2s + z)C(« - z) 


(3/2) 

C(*)C(-* 


r(s) 


dz 


dz = - 
zsiti{ttz) 2 


(3/2) 


C (| + #) C (~| -it) 

(| + it) cosh(-7rf) 


dt, (57) 


since sin (ir (| + if)) simplifies to — cosh(Trf). Combining (57) with the com¬ 
putations leading up to (56) hnishes the proof of Theorem 1.4 □ 


6 Asymptotics of generating functions for SU( 3) 
representations 

In this section we begin applying our results on the function ur(s) to derive 
asymptotic results for the enumeration of representations of SU( 3). Let us 
start with a few general remarks about the broader context for our analysis. 
Theorem o is formulated most naturally as an asymptotic enumeration 
result for representations of SU( 3), but one can think of the sequence r(n) 
as enumerating a certain combinatorial class of integer partitions, namely 
the class of representations of an integer n as a sum of integers of the form 
Qj,k = \ jk{j + k), where repetitions are allowed; the order of the summands 
is unimportant; and if = ay ^ but (j, k) / ( j',k'), the numbers ajy 
and ay k' are considered distinct. Thus, if one ignores the representation- 
theoretic aspect, our result broadly falls into the area of partition asymp¬ 
totics. This topic, which had its beginnings in 1918 with the asymptotic 
expansion of Hardy and Ramanujan [20] for the partition-counting function 
p(n), has important connections to analytic number theory and the theory 
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of modular forms (see ED- Highlights of the theory include Rademacher’s 
refinement of Hardy and Ramanujan’s formula into a convergent series for 
p(n) j40| and Wright’s asymptotic formula |49| for the number of plane par¬ 
titions of n. See 0, m, EE 121 for more recent results. 

One landmark result in the theory was the theorem of Meinardus, who 
in his 1954 paper m identified a general strategy for proving asymptotic 
formulas for a wide class of combinatorial enumeration sequences c(n) whose 
generating functions have an Euler product formula of the form 

OO oo 1 

1 + £ c(n)x n = [] (1 _ xm)am (58) 

n=l m= 1 ' ' 

for some sequence (ct m )m=i °f nonnegative integers. Assuming certain tech¬ 
nical conditions are satisfied, Meinardus’s result states that c(n) has an 
asymptotic formula of the form 

c(n) = (1 + c^l^CVr -7 exp (59) 


for certain explicitly computable constants C, D, 7 , 6. His result includes 
many earlier results (including ([Tj) and Wright’s formula for plane partitions 
mentioned above) as special cases that can be deduced from it without 
much difficulty; see [31 Chapter 6 ] for examples and an accessible exposition 
of Meinardus’s result. 

Now note that the generating function ([ 2 ]) of our sequence r(n ) can in 
fact be written in the form (58), with the exponents a m being given by 
a m = #{(j, k ) £ N 2 : m = jk(j + k)/2} (that is, the number of inequivalent 
m-dimensional irreducible representations of SU(3)). One might hope that 
Meinardus’s theorem will apply in this case, but it turns out that it does 
not, and indeed, the asymptotic formula ([3]) ends up having a more elaborate 
form than (59). Nonetheless, the general strategy employed by Meinardus 
and earlier authors for proving asymptotic formulas such as 0 and ( |59| ) 
(described for example in |15( Section VIII. 6 ]) remains valid, and we were 
able to adapt it to our needs, although a few nontrivial technical hurdles 
need to be overcome. One conceptual innovation that simplifies the anal¬ 
ysis somewhat is to use a probabilistic representation similar to Fristedt’s 
probabilistic model for random integer partitions m- 

The first step in the analysis, which we undertake in this section, is to 
understand the asymptotic behavior of the generating function of r(n) and a 
few related functions, as the argument approaches a singularity point. The 
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main tool we will use is the Mellin transform and its inverse. Define 


oo _. oo 


G{x) 

n i _ xd fc (i+ fc )/2 

r(n)x n 

(1*1 

< 

1) 


i fc=i 

n= 0 




fit) 

oo 

= exp ( - lyT’(y + 

A:)f) 

if 

> 

0) 







h(t) 

= log G(e _t ), 


it 

> 

0) 

Note that the 

trivial bound 






( 60 ) 

(61) 

(62) 


OO oo 

x yk{j+k)/2 < ^ | x |(j+A;)/2 

j,k =1 j,k= 1 




2 


shows that the product defining G(x) (hence also the series) is absolutely 
convergent for complex x satisfying \x\ < 1 and defines an analytic function, 
and that /(t), h(t) are defined and finite for t > 0. We will also need a few 
additional easy bounds that are given in the following lemma. 

Lemma 6.1. The functions f(t) and h(t ) satisfy the asymptotic bounds 


m 


o(r 2 / 3 ) as t \ 0, 
0{e~ t ) as t —»• oo, 


h{t) 


0(t~ 2 / 3 ) as t \ 0, 
0{e~ t ) as t —> oo. 


Proof. Consider first /(f); the bound /(f) = 0(e~ t ) as t —»• oo is trivial and 
left as an exercise. The behavior of /(f) as t \ 0 can be understood by 
writing 


m = r 2 / 3 V e xp (—J(i 1/3 j)(i 1/3 fc)(i 1/3 j +t 1 / 3 *) 4 123 1 '/ 3 , 


and noting that this is f -2 / 3 times a Riemann sum (with Ax = Ay = f 1 / 3 ) 
for the double integral I = exp ( — \xy(x + y)) dxdy. This integral 

can be seen to be finite by making the two-dimensional change of variables 
u = x, v = \xy[x + y), which, after a short computation that we omit, gives 


n OO | 

e~ v = dudv. 

Vu 4 + 8 vu 
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Making another substitution, namely w = ( 8v ) 1 ^ 3 u for the integral with 
respect to u, transforms this into 



proving therefore that /(f) = Oft 2 / 3 ) (in fact f(t ) = (I + o(l))t 2 / 3 ) as 
f \ 0. Finally, for h(t), note that 


h(t) 


- lo g ( x ~ ex p + *0*) ) 

OO OO ^ 

EE — exp + k)mt) 

j,k =1 m=l 


CAJ I » CAJ 

= s - 5 Z ex p (-yw+ k ^ mt ) = n -/m)- ( 6s ) 

m= 1 \j,k= 1 j m= 1 

The claim about the behavior h(t) as t —> oo is again trivial, and for t near 
0 we have 

°o / 00 1 \ 

|h(f)| < ^ —l/M)l = O ^ —(mf)“ 2/3 

^ m \ ' m I 

m =1 \m=l / 

= 0 (E ^73 ■ t_2/3 ) = °< r2/3 >' 

\m=l / 

(Note that the 0(t~ 2 / 3 ) bound for /(f) actually holds for all t > 0, since e _t 
decays faster than t -2 / 3 as t —> oo.) □ 


We are now ready to prove a much stronger result about the asymptotic 
behavior of h(t) and its derivatives as t \ 0. 

Theorem 6.2. As t \ 0, the function h(t ) has the asymptotic expansion 
h(t ) = pi_t 2 / 3 + fj, 2 f -1 / 2 _ i logf + ca'(0) + | log 2 + O^ 1 / 2 ), (64) 

where m,p 2 are constants defined by 


_ ^'p/n 2 , 5 , 

V 1 — 3 1 U) £ U) 
h2 = V^C (2) C ( 1 ) • 
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Furthermore, the expansion (64) can be differentiated termwise; specifically, 
the first two derivatives h!{t) and h"(t) have the asymptotic expansions 


h!(t) = — ^nit 5//3 — \ 1 i 2 t 3//2 - \t 1 + 0(t 1 Z 2 ), (65) 

h"(t) = f Mii" 8/3 + l/xa ^ 5/2 + £i" 2 + 0(t- 3/2 ). (66) 

Proof. First, we wish to show that f(t), and therefore also h(t), can be re¬ 
lated to the analytic function u>(s) via the Mellin transform. Start by consid¬ 
ering f(t). Computing its Mellin transform, which we denote by A4[f(t)](s), 
we have that 


/»oo txj poo 

M[f{t)](s) := / f(t)t s ~ 1 dt= E / exp (- \jk{j + k)t)t s ~ l dt 
Jo j.k=i Jo 

OO 

= E (iW + * 0 ) _S r ( s ) = 2 T(s)w(s). ( 67 ) 


Next, consider h(t). Taking the Mellin transform of (63), we obtain easily, 
using the standard linearity and scaling properties of the Mellin transforms, 
that 


OO 1 OO ^ 

M[h(t)](s) = E -M[f(mt)fis) = £ -m~ s M[fms 

z ' m Z —' m 


m =1 


m= 1 


= 2 s r(s)cj(s)C(s + 1 ). 


( 68 ) 


By Lemma 6.1, both (67) and (68) are valid for all complex s satisfying 
Re(s) > 2/3. (This is important for the proof below but also provides a 
useful consistency check with Theorem mi)-) 

We are now in a good position to apply the method of Mellin transform 
asymptotics, described for example in m Appendix B.7]. The main idea is 
that the asymptotic behavior of h(t) as t \ 0 is closely tied to the singularity 
structure in the complex plane of its Mellin transform. Specifically, we apply 
the Mellin inversion formula ESI eq. (3.1.5), p. 80] to deduce that h(t) has 
the contour integral representation 


h(t) = 


1 

2^xi 


2 s T(s)o;(s)C(s + l)t s ds, 


(69) 


(“) 


where a > 2/3, and (a) denotes as before the vertical contour from a—ioo to 
a + ioo. Shifting the contour of integration to the left past a few of the poles 
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of the integrand will now produce the main asymptotic contributions from 
the residues, and the integral will become an error term whose magnitude 
is easily controlled. The relevant poles past which we will move the contour 
are at s = 2/3,1/2,0, —1/2. The first two are simple poles, with respective 
residues 


RgS s= 2/3 

Res s=1 /2 


2 s T(s)o;(s)C(s + 1 )t 


2 2/ 3 


r(l) 2 C (|)t 


5\ +—s 


2 s r(s)w(s)C(s + l)t s = %/2/rC (I) C (|) t 


Hit s , 

H2t~ s - 


The pole at s = 0 is a double pole since each of the factors T(s) and ((s + 1) 
contributes a first-order singularity. By the well-known facts that T(s) = 
1/s — 7 + O(s) and C(s + 1) = 1/s + 7 + O(s) near s = 0 (where 7 is the 
Euler-Mascheroni constant), the Laurent expansion of 2 s T(s)o;(s)C(s + l)t~ s 
around s = 0 is easily computed as 


T(s)o;(s)C(s + l)(t/2) s = - - 7 + O(s) - + 7 + 0(s) 


x (i+ u}'(0)s + 0(s 2 )) (1 - s log(t/2) + 0(s 2 )) 

= ^2 + ( w '(°) + l lQ g 2 - l l °St) \ + 0(1). 

Thus, the residue at s = 0 is u/(0) + ^ log 2 — | logt. Finally, the residue at 
s = -1/2 is (jy / 7r72C(-5/2)C(l/2))t 1 /2 =: 

Now shift the contour in ([69]) past the four poles, say to ( a') where 


a' = —1. By the residue theorem (with easy justification provided by (45), 
@ and ([it])), we get that 


h(t) = nit 2 0 + Hit 1//2 — | logt + ^(O) + | log 2 + ut 1 ^ 2 

+ ——: f 2 A T(s)o;(s)C(s + l)f s ds. 

2m J 

(- 1 ) 


(70) 


It is easy to see that the integral is 0(t 1 ), so this proves (64). 
{65)—(66) follow by differentiating both sides of (70). 


Finally, 

□ 


7 A probabilistic model for representations of SU (3) 

Next, we take advantage of the product structure of the generating function 
© to interpret the coefficients r(n) in terms of a probabilistic model involv¬ 
ing a family of independent random variables distributed according to the 
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geometric distribution. Such probabilistic models are widely applied to the 
asymptotic enumeration of combinatorial structures and to the analysis of 
their probabilistic properties (see ED- A version specifically tailored to the 
study of integer partitions was introduced by Fristedt [16] (and is sometimes 
referred to as Fristedt’s conditioning device) and forms the prototype for a 
variety of similar models, considered, e.g., in [45] • 

Let t > 0 denote a parameter. For integers j. k > 1, let Xj k denote a 
random variable defined on some probability space having the distribution 
X ~ Geomo(l — e _:,fc ^ +fc ^ //2 ), namely the geometric distribution (in the 

version that starts at 0, hence the notation Geomo(-)) with parameter P j k = 
1 _ e -jk(j+k)t/2 _ That ig) 

we have 

p t(X jtk = m) = = (1 - e -Mj+k)tl2 )e -M^)mt/2 

for m > 0, where the notation Ft(-) denotes the probability with respect 
to the parameter value t. We assume that all the variables (Xj k )j k >i are 
simultaneously defined on the same probability space and form an indepen¬ 
dent family of random variables. 

We interpret the variables (Xj }k )j tk >o as encoding the structure of a 
certain random representation of SU( 3), namely 


W = 


0 


j,k 


(71) 


j,k =1 


which consists of a sum of Xj k copies of each of the irreducible representa¬ 
tions Wj' k (see Subsection 1.1). The total dimension of W is therefore the 
random variable, which we denote N, given by 


N = 


X! \iKj + k)X j)k . 


(72) 


j,k= 1 


It is easily seen from the Borel-Cantelli lemma from probability theory that 
with probability 1 only finitely many terms in the sum are nonzero, and 
therefore that N is almost surely finite. 

If {rnj )k )X k= i is an array of (nonrandom) nonnegative integers, we have 


n x* 

VL fc=1 

„—nt 


G{er*Y 


m j,k} 


= _ e -jk(j+k)t/2 

j,k =1 


e -jk(j+k)mj, k t/2 


(73) 
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where G(x) is defined in ( f60] ) , and n = ^°° fc=1 \jkij + k)m,j } k■ This can be 
interpreted as the probability in our random model of the event that the 
random representation W in ( |71[ ) is equal to the specific nonrandom repre¬ 
sentation ®’j°k=i rn j, k ^j,k^ which has dimension n. Note that representations 
with the same dimension have equal probabilities of being observed under 


the measure IPt(-). Summing the probabilities (73) over all r(n) representa¬ 
tions of dimension n therefore gives that 


= n) = 


e nt r(n ) 


G{e 




(74) 


This key formula will be of crucial significance in our analysis; it relates 
our integer sequence r(n) to the distribution of the random variable N, 
and will therefore allow us to translate results about that distribution into 
information about r(n). 

Let IEt (•) and Var t (-) denote the expectation and variance, respectively, 
relative to the parameter value t. The following lemma provides another 
connection of the random variable N to a function we have been studying. 


Lemma 7.1. The expectation and variance of N are given, respectively, by 


Ei(iV) = -h\t), 

\ait(N) = h"(t). 

Proof. Denote g{t) = G{e ~ t ) = 0 r(n)e~ nt , so that 

OO OO 

g'(t) = -^2 nr(n)e~ nt , g"(t) = ^ n 2 r{n)e~ nt . 

n =0 n =0 


Using (74), it follows that 


h'(t) = (log g)'{t) = 


g'if) 

gif) 


Y.nPtiN = n) = -Et(N), 

n =0 


and similarly, 

g"{t)g(t)-g\t? _ g"(t) (gft)\ 2 

git ) 2 git) 

OO 

= J2 n 2 Wt(N = n)~ i-E t (N)) 2 = E t (N 2 ) - (E t (N)) 2 = Var t (JV). □ 

n =0 
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8 Saddle point analysis 


Our next step is to perform a saddle point analysis of the generating function 
G(x). This can be described in two equivalent languages. One standard 
description of the method commonly found in analysis textbooks is in terms 
of the complex-analytic idea that the coefficients r(n) may be extracted from 
G(x) as contour integrals, namely as 


r(n) = - 

v ' 2m 


G(z ) 

?n+l 


dz. 


(75) 


\z\=r n 


The idea is then to choose the radius r n of the contour of integration in 
such a way that the bulk of the contribution to the integral comes from the 
environment of the point z = r n on the positive x-axis, and the nature of 
this contribution can be understood from the local asymptotic behavior of 
G(x) near x = 1. 

An equivalent way to describe the saddle point technique (as applied to 
the present context) is in probabilistic terms, based on the relation (74). The 


idea is that we are free to choose the value of the parameter t in the random 
representation model, and we will do so in a way that causes the random 
variable N to have its mean value at (or near) n. Probabilistic intuition 
then predicts that N, being the sum of independent random variables, will 
have an approximately Gaussian distribution, and from this intuition an 
asymptotic formula for the left-hand side of © (hence for r(n)) can be 
easily guessed. This guess, known as a “local central limit theorem,” can 
then be proved using Fourier inversion, which is formally equivalent to the 


contour integral (75). Despite the equivalence of the two approaches from 


the point of view of mathematical analysis, they are not psychologically 
equivalent, and we prefer the probabilistic approach for putting the analysis 
on a more conceptual footing. 

Say that a sequence of positive numbers (t n )^ = i is an (approximate) sad¬ 
dle point sequence for the random representation model if the asymptotic 
relation 

E tn (AT) = -h\t n ) = n + 0(n 7 / 10 ) (76) 

holds as n —> oo. (As often happens with saddle point analysis, it is not 
crucial to identify the saddle point precisely, and for reasons that will be¬ 
come apparent shortly, an approximate expression with 0(n 7 / 10 ) error will 


be sufficient for our needs.) Our first goal is to find solutions to (76). Re¬ 


ferring to the asymptotic expansion (65), we see that in order for — h'{t n ) to 


grow to infinity, as it must according to (76), the sequence t n will decrease 
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towards 0. Considering only the leading term in (65) temporarily, we get an 
approximate equation 


|/iit n 5/3 


n. 


which can be easily solved to show that t n ~ (3n/2//i)~ 3;/5 . That is, t n 
should decay roughly as a constant times n _3//5 , where the constant is 
(|/ri ) 3 / 5 = 2X 2 (with X as defined in Q). As it turns out, this expression 
is not sufficiently precise for our needs because the contribution from the 
second-order term in (65) (proportional to t n 3 ^ 2 x n 9//10 ) is too large, 


so we 


look for a more precise solution to (76) whose form is given by the ansatz 


in = rin - 3 / 5 -T 2 n- 7 / 10 -T 3 n- 4 / 5 , 


(77) 


with n = 2X 2 and t 2 ,T 3 being constant coefficients whose value needs to 
be determined. Solving for r 2 and 73 is now an amusing calculus exercise. 
Introduce symbols q, u, v, defined by 


q = 


n -1 / 10 , 


Then t n can be written as 


T2 

T3 


V = — 

n 

n 


t n = ng 6 ( 1 -uq- vq 2 ) = rig 6 (l - q(u + vq)). 


(78) 


(79) 


Now consider each of the first two leading terms in ( |65[ ) with t = t n . Using 
the second-order Taylor expansion (1 — z) a = 1 — az + \a{a — 1 )z 2 + 0(z 8 ) 
and expanding in powers of q , we have 


im * n 5/3 = ( r i^ 6 ) / {} ~q(. u + vq)j 

= «- 10 (l - (- 1 ) ?(« + vq) + \ (- 1 ) (- 1 ) q\u + vq? + 0{q?) 
= q~ 10 + (fit) q~ 9 + (§« + §■ | u 2 ) q~ 8 + 0(q~ 7 ), (80) 


and similarly, 

lfJ-2tn 3/2 = 1^2 ( T l<? 6 ) 1 (i - q( u + v q )) 

= (i/i 2 u“ 3/2 ) 9 _9 (l - (-|) 9(« + v?) + 5 (- 1 ) (-|)^ 2 (« + ^) 2 ) 
= (^ 2 Ti _3/2 ) q~ 9 + (| • l/J, 2 T? 3 / 2 u) q~ 8 + 0{q~ 7 ). (81) 
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The remaining terms are 0(t n 1 ) = 0(q 6 ) = 0(n 6//1 °), which is within our 
tolerance range for solving (76). Thus, adding up ( [80] ) and (81), we get that 

=n + (§ u + n 9/10 


+ ii«+i-i « 2 


+ 1 ■\^rf /2 u) n 8 / 10 + O(n 7 / 10 ). 


Comparing this expansion to (76) shows that the sequence (t n )™ =1 will be a 
saddle point sequence provided that the equations 


5 1 -3/2 

3« = -2/i2r 1 , 

5 4 2 3 1 


V ~ 3 ' 3 U 2 2^ 2T l 


--3/2 


U, 


hold, or equivalently, provided that T 2 ,r 3 take the values 

3 — 1/2 
72 = -100271 , 

4 2-1 9 -3/2 

T 3 = -gTaTi - 20/X2T! T 2 . 

After a small amount of further algebraic simplification, which we omit, 
we arrive at the following result. 

Lemma 8.1. Let the numbers X , Y be defined by © ©■ Define a sequence 
of numbers {t n )^ =1 by 

t n = Tin -3 / 5 — T 2 n^ 7 / 10 — T 3 n~ 4/l5 (n > 1), (82) 


where n, r 2 , 73 are the constants 


n = 2A , 


T9 = —-Y _1 y, 

10 

r 3 = —A- 4 T 2 . 
400 


Then (t n )^_ 1 is a saddle point sequence, that is, it satisfies the equation 


(76). 


It is worth noting that the expression (82) for the location of the saddle 


point requires considerably more precision than the case of Meinardus’s 
theorem on partition asymptotics (discussed in Section [6]). In Meinardus’s 
analysis the saddle point could be defined using a single term, given as a 
constant times a power of n, analogous to the leading term nn -3//5 in (82); 
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see equation (6.2.15) on page 93 of jlj. The reason is that Meinardus’s 
theorem makes the relatively restrictive assumption that the Dirichlet series 
Ylm a m'm-~ s (in the notation of (58)) can be analytically continued to the 
strip Re(s) > <7o for some gq < 0, with only a single pole of order 1 at some 
point s = a > 0. In our situation the relevant Dirichlet series is 2 s w(s), 
which has not one but two separate poles with positive real parts, namely 
at s = 2/3 and s = 1/2. It is precisely this somewhat novel singularity 
structure that leads to the more involved saddle point analysis leading up to 
(82), and which (as we shall see in the next few sections) ultimately causes 
the asymptotic formula ([3]) for r(n) to have the interesting structure it does. 


9 Saddle point analysis II: the local central limit 
theorem 


Lemma 9.1. As n —oo, the variance of N at the parameter value t = t r 
is given asymptotically by 


Var tn (N) = (§X- 2 ) n 8 / 5 + 0(n 3 / 2 ). 


(83) 


Proof. Since, by Lemma 7.1, the variance is given by Var t n {N) = h"(t n ), 
substituting the values (82) into the asymptotic expansion ( |66| ) gives the 
result after a quick calculation, which we omit. □ 

We are now ready to formulate our local central limit theorem for the 
random variable N , whose proof will be the heart of the saddle point anal¬ 
ysis. Note that the saddle point equation (76) identifies the expe cted value 

shows 


9.1 


E t n (N) to within an error of 0(n 7//1 °). On the other hand, Lemma 
that N has standard deviation (relative to the parameter value t = t n ) of 
order n 4 / 5 = n 8 / 10 . The local central limit theorem is the statement that 
the distribution of N is asymptotically Gaussian near its mean value—or 
equivalently, near the value n, since the two values are 0(n -1 / 10 ) standard 
deviations apart (this was the reason we required an error no greater than 
0(n 7 / 10 ) in (76)). The precise statement is as follows. 


Theorem 9.2 (Local central limit theorem for N ). The variable N “satisfies 
a local central limit theorem near its mean value. ” More precisely, as n 
we have the asymptotic relation 


oo 


\(N = n) = 


l + o(l) 


V^vrVar^AO 


= (1+ »(D) 

V 5vr 


(84) 
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It is well-known that local central limit theorems for integer-valued ran¬ 
dom variables are inherently a more delicate and less robust phenomenon 
than the usual, non-local sort, and this case is no exception. The point is 
that the definition of N involves a sum of independent, but not identically 
distributed, components \jk{j + k)Xj^. each of which takes its values in 
the sub-lattice ( \jk{j + k))Z of Z. The proof of local CLTs of this type 
therefore invariably requires ruling out the potentially harmful influence of 
periodicities. (For example, in the most trivial example illustrating this 
principle, if the random variable we were looking at were a sum of random 
variables taking even values, then there would be zero probability for it to 
take on an odd value; however, even the absence of rigid or “determinis¬ 
tic” periodicities of this sort leaves room for more fuzzy periodicities of a 
probabilistic kind.) In our setting, this is at heart a number-theoretic claim, 
and the techniques of analytic number theory are the most appropriate to 
use to attack it. Specifically, it turns out that the modular transformation 
properties of the Jacobi theta functions are sufficient to prove the bounds 
we will need. 

Before proceeding to read the proof of Theorem 9.2 which we present 
in this section and the next one, the reader may wish to skip directly to 
Section elu where we show how Theorem |9.2| and the other results proved 
so far imply our theorem (Theorem 1.1) on the asymptotic enumeration of 
representations of 517(3). 


The proof of Theorem |9.2| begins with the characteristic function (a.k.a. 
Fourier-Stieltjes transform) of N, which we denote 4>t{u) and is given for 
u e M by 


Mu) = E t (e iuN ) = ^>t(7V = n)e inu 

n =0 


OO / \ —Tit 

V- r(n)e nt inu 
M G(e~ t ) 

r> —D ' ' 


1 

W *) 


oo 

£ 

71=0 


r{n)e~ nt+inu 


G(e ~ t+iu ) 
G(e-*) 


(85) 


(As before, the subscript t emphasizes the dependence on the parameter t, 
which soon we will specialize to the saddle point values t = t n .) Since N is 
an integer-valued random variable, this is simply a Fourier series. Fourier 
inversion therefore gives that 

P t (N = n) = M Mu)e~ inu du. ( 86 ) 


Denote a n 


(Vari ?1 (-/V)) 1 / 2 , the standard deviation of N. Setting t = t n 
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and scaling the integration variable u by l/a n in (86), we get 


Ptn (N = n) = 


2-7TOV 


(fi tn (u/a n )e~ inu ^ du. 


Comparing this to the claim (84), we see that it will be enough to prove 
that 


^t n {u/a n )e mu/rTn du 


V2tt. 


(87) 


A key first step will be to prove pointwise convergence of the integrand 
to the characteristic function of the standard normal distribution. 


Theorem 9.3. For each «£R, we have 

<t>t n {u/Gn)e - inu - 


->• e 


—v?/2 


( 88 ) 


Note that Theorem |9.3| is equivalent to the claim that the normalized 
random variable N = a~ 1 (N — Kf n (N)) converges in distribution to a stan¬ 
dard normal random variable, that is, to a non-local CLT. 

The proof of Theorem 9.3 involves not much more than a Taylor expan¬ 
sion, but one needs to make sure that the dependence on both the parameters 
t and u is carefully taken into account. We precede the main part of the 
proof with two technical lemmas. To motivate the statement of the first 
lemma, recall that if A is a random variable with the geometric distribution 
Geomo(p), where the parameter 0 < p < 1 is expressed as p = 1 — e _t , then 
the mean and variance of X are given by 


E(X) = 


1 ~P 
P 


t 


1 — e 


-t ■ 


Var(X) = 

pZ 




(1 — e - *) 2 ’ 


and the characteristic function (j)x(u) = E(e*“ A ) of X can be easily com¬ 
puted to be 


00 i —t 

fa(») = £u - e-'K'V” = . 

n =0 

In particular, if log z denotes the principal branch of the logarithm function, 
then log (fix (it) has the Taylor expansion 

log i fix(u) = i E(A)n — | Var(A)u 2 + 0(u 3 ) as u —> 0. 
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Lemma 9.4. Define a function 




—t 


T 1 {fiu)=\0g l l _^ t+iu ~i 

For some constant C > 0, the bound 


i 37« + 5 

1 — e * 




(1 -e 


—t\2 


U 


(t > 0, it E 


\v(t,u)\ < C 


e t \u\ 3 
(1 - e _t ) 3 


holds for all u £ M and t > 0. 


Lemma 9.4 has been previously used in the study of integer partition 
asymptotics. See pp. 10-11 of [12] for the proof. 

Lemma 9.5. The improper two-dimensional integral 

roc roo / xy ( x + y jS\3 exp {-\xy(x + y)) 


I = 


>0 Jo 


{1 - exp (~\xy{x + y)))' 


dx dy (89) 


converges. 


Proof. Using the expansion = Yl^= l to expand the integrand, 

we get that 




n + 1 


n= 1 
oo o 


/O 40 


+ y) 


exp (— ^nxy(x + y)) dx dy 


< 


v /*oo /»c 

L 


ab{a + 6) 


exp (— ^ab(a + 6 )) dadb 


/*oo /*c 

= C(|) / 

Jo Jo 


ab{a + b ) x 3 


exp (—|a 6 (a + 6 )) dadb. 


The fact that this last integral converges can be shown easily using the same 
change of variables used in the proof of Lemma 6.1 D 


Proof of Theorem 9.3. Using (85) and (60), we represent 4>t n (u/a n )e inu ! an 
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as an infinite product. Taking the logarithm, we have that 


log 


h n (u/a n )e- inu ^ 


= log 


G(e *" + * u / <T ") inu 


G{e 


tn 


(Jr 


1 - exp ( -\jk{j + k)t n ) 


b \ 1 - exp {\jk(j + k)(-t n + iu/a n )) 


mu 


E 

j,fc=l L 


jk{j + k)t n jk(j + fc)u 
2 ’ 2 a n 


+i 


. \jk(j + k) exp {—\jk{j + k)t n ) u 


1 - exp (■ -\jk(j + k)t n ) 


(Jr 


i (\jKj + fc)) 2 exp (~|jfc(j + k)t n ) ur 
(l-exp (~\jk{j + k)t n )) 2 a n 


mu 


(J 77 


(jk{j + k)t n jk{j + k)u\ . 






— ^ Vai't„(iV)^ = — Tjtt 2 + 0(n 1//10 u) + i? n (n), (90) 


(T“ 


where we denote R n (u ) = ^°° fc=1 V y 2 
this quantity can be bounded as 


( jk(j+k)t n jk(j+k)u 
’ 2<r n 


Invoking Lemma 


9.4 


\Rn{u)\ < ^ 

j,k=1 


/ jk(j + fc)t ra jk(j + k ) u 


< c l - 


\ 2 

jk{j + A:) 


2(7" 71 


r 


■E 

j,k =1 
, 13 00 


exp {-\jKj + fe)^n) 

(l - exp ( ~\jk{j + k)t n )Y 


= C 


3,11/3 

j,k =1 


Wfc(j + fc)\ 3 exp {—\jk{j + fc)f n ) . t 2 /3 
(l - exp (~\jk(j + k)t n )) 3 


In this expression, the prefactor Clnl 3 ^ 3 ^ 11 ^ 3 is of order 0(n 1 / 5 |n| 3 ), 

1 /3 

and the sum is a Riemann sum (with Ax = Ay = £ n ' ) for the double 


integral (89). Since that integral is finite by Lemma 9.5, the Riemann sum 


converges to it and in particular is bounded. This together with the above 
computations establishes that 


log 


$ tn (u/<r n )e- inu ^ 


-> -in 2 , 
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which implies ( 88 ) and therefore finishes the proof. 


□ 


The computations in the proof above also imply a useful bound on 
\<t>t n (u/a n )l which we summarize in the following lemma. 

Lemma 9.6. There exists a constant D > 0 such that the bound 

Wt n {u/a n )\<e- u2 ' A (91) 

holds for u satisfying |u| < Dn 1 ^ 5 . 

Proof. By our earlier computations, we have that 


l&r.faMOl = exp Re^log 4>t n (u/ cr n )^j < exp (^-\u 2 + Cin 1/5 |u| 3 ) 


for some constant C\ > 0. (Note that the term i{K tn (N) — n)u/a n in (90) is 


a pure imaginary number and therefore only contributes a phase factor to 
4>t„(u/a n ).) Restricting u to satisfy |u| < Dn 1//5 for some sufficiently small 
D > 0 guarantees that — \u 2 + C'in~ 1 /' 5 |u | 3 < — \u 2 , so in that range we get 
the bound ((Ml). □ 


10 Saddle point analysis III: the Jacobi theta func¬ 
tions 


In order to deduce (87) (and hence Theorem 9.2) from Theorem 9.3 we will 
need the “easy” tail bound ( pH] ) together with a much more delicate bound 
on \4>t n (u/cr n )\ that can be shown to hold outside the range \u\ < Dn 1 ^. 
The goal of this section is to prove this more difficult bound, given in the 
following theorem. 


Theorem 10.1. There exists a constant E > 0 such that the bound 

\(ft n (u/a n )\ < exp (-£n 3 / 10 ) (92) 

holds for u satisfying Dn 1 ^ < |w| < ira n , where D is the constant from 
Lemma 1.9. (A 


As preparation for the proof, replace u/a n by u, so in view of the asymp¬ 
totics (83) of the variance, the statement we need to prove becomes that 


\f>t n {u)\ < exp (-EA 3 / 10 ) 


(93) 
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for some constant E > 0, all sufficiently large n, and all u such that 
D'n ~ 3 / 5 < |tt| < 7 r, where D' = 21A/(y / 5/6X~ 1 ). More generally, we will 
prove that for any e > 0, there exist numbers (3,6 > 0 such that if 0 < t < (3 
and et < lul < 7r then 


\M U )\ < exp (~5t 1/2 ) 


(94) 


The bound (93) then follows by substituting t = t n (refer to (82)). 
Now observe that 


logcj)t(u) = logG(e t+ *“) — log G(e f ) = h(t — iu) — h(t) 

V . 

= X] — (/M* - iu )) - /M))- 

z —' m 

m=l 


By the definition of f(t) in (61), clearly R e(f(t 
u £ M, t > 0, so 


iu)) < R e(f(t)) for any 


Re( log^t(u)) < ^2 — Re(/(m(f - iu)) - f(mt)) < R e(f(t - iu) - f(t)) 

7Tb 

m =1 


< Re 


^ (exp + k)(t - iu)) - exp + k)t) ) 

j,k =1 


In this last expression, the real part of each of the summands is nonpositive, 
so we can omit as many of the terms as we like and still get an upper bound. 
We focus on the subset j = 1 of the summation range, leading to the bound 


Re (log Mu) ) < Re 


exp [—^k(k + 1 )(t — iu)) — exp (y—\k{k + 1 )t) 


k =l 


k =1 


= Re(A(i — iu)) — A (t), 


(95) 


where we denote 

OO 

A(z) = ^ exp (— \k(k + l)z) . 
fe =0 

Proposition 10.2. For any e > 0, there exist numbers (3, 7 > 0 and 
0 < p < 1 such that for all t,u satisfying 0 < t < (3 and et < |«| < it, 
the bounds 

m > 7 t- i/2 , (96) 

hold. 
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Proposition 10.2 is enough to prove (94), and hence Theorem 10.1 


since 


(95) and (96) imply (under the assumed conditions on t,u) that 

\M U )\ = exp (Re(log H u ))) < exp ( - A (t) + (1 - p)X(t)) 

= exp(— pX(t)) < exp . 

To prove the proposition, we will use the fact that A (z) can be related 
to the Jacobi theta functions. Define 


(k+ 1/2) 2 = ,„l/4 ^ 0 k{k+ 1) 
k =0 
k‘ 2 


02 (q)= J2 q {k+1/2)2 =2q l/A Y J ( l 

k=—oo 
oo 

Hq) = q 

k=—oo 
oo 

«4(4)= Y, MJVA 
k =—oo 


H*) = e 2 (e™ z ), 
Hz) = He inz ), 

Hz) = He inz )- 


The functions J 2 (z), Js(z), Ja{z) are holomorphic functions in the upper half 
plane H = {x + iy : y > 0}. They are known to satisfy the modular 
transformation relations 


Hz) = (- iz ) 1/2 H~ 1 /z), (97) 

H z ) = (-**r 1/2 J 3 (-V*), (98) 


where w H > w 1//2 is the principal branch of the square root function; see m 
p. 181]. In particular, note that for real y > 0, J^(iy) = 1 + 0{e~ ny ) and 
Ji(iy) = 1 + 0(e~ 7Ty ) as y —> oo. Therefore (97) and (98) imply that as 
j/-l0we have 


Hw) = Yj h{ - i/y) = ^ ( x + 0 ( e_7r/ ")) ’ (") 

Hw) = ^Hi/y) = ^ (i + o (e^)). (100) 

Furthermore, we also have trivially that | J^{x + iy)\ < Jz{iy)- 

Now we can relate the previous discussion regarding A (t — iu) to the 
functions J 2 (z), J-s(z), Ji( z )-> by noting that 


OO 

A (t — iu) = ex P + l)(t — iu)) 

k =0 


l e t/S e -iu/Sj 2 


u + it A 
2 vr ) ' 
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In particular, when f \ 0 we have the bound A(t) = J2(it/2ir) = 
^ ^1 + O y / 27r/t > 7 i -1 / 2 , which was the first claim of Proposi¬ 

tion 10.2 For the second claim, by modifying the constants we can forget 
about the 1/tt factor, and prove instead that 

- iu))\ = \J2 (^)| 

A(- 7 r t) J2(it/2) 

is bounded away from 1 (more precisely, is < p for some number p < 1 ), 
under the assumption that: 

1 . et < | it | < 1 where e > 0 is given; 

2 . 0 < t < f3 where we are free to fix /3 as small as we please. 

(Naturally, j3 and p will depend on e.) Observe therefore that, thanks to 
we have 


MU/2) l(„4.«W9l 1 /2 l T V ^ 


J 2 {it/ 2 ) \[u + it)/2\ 

< (l + 0(e- 27r /*)) ^ 


J 4 - 


2 u 


u* +1 2 


+ i- 


21 


+ t 2 


1/4 


U‘ 


+ t 2 . 


■h ( * 


2 t 


+ t 2 


U‘ 


( 101 ) 


To bound this expression away from 1, we divide into three cases: 


Case 1. Assume that et < |it| < t 3 / 4 . In this case, 


t 2 + u 2 


1/4 


< 


if 


t 2 + e 2 t. 2 


1/4 


1 + e 2 


1/4 


The quantity on the right-hand side is smaller than 1, and we denote it 
1 — p\. Second, making the further assumption that t < 1 (which simply 
puts the constraint /? < 1 on our ultimate choice of /3), we have 

2t . 2t . 2t 1 


> 


> 


U 2 + f 2 “ t 3 / 2 + t 2 “ 2f 3 / 2 V*’ 

so that (by the bound for Jz(iy) as y —>• oo) J 3 727 / 7 ^ = 1 + 0{e~' Kt i/2 ) 
as t —>• 0. Combining the last two estimates, we get that the right-hand side 
of ( 101 ) is bounded by 

(l + O (e- 2 ^)) (l + O (e"^)) x (1 - Pl ), 
which is bounded by 1 — pi /2 if /3 is assumed to be small enough. 
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Case 2. Assume that f 3 / 4 < |u| < Mf 1 / 2 , where M is some large enough 
constant whose value will be specified shortly. In this case, note that, again 
using the assumption that t < 1, we have 

2t ^ 2t 21 _ 2 

u 2 + t 2 ~ M 2 t + t 2 ~ (M 2 + l)t “ M 2 + 1' 

The function s t —> Js(is) of a positive real variable s is monotone decreasing. 
So, after fixing M , we get that J 3 is bounded by L := J 3 

in the range of values under discussion. Furthermore, we have 


f 2 + u 2 


1/4 


< 


t 2 +f 3 / 2 




1/4 


t l/8 


By restricting (5 (and therefore t) to be less than 



we ensure that 


t 2 

t 2 + u 2 




< t 1/s L < 1/2, 


which, when combined with (101), proves the claim in this case. 


Case 3. Finally, assume that Mt 1//2 < |u| < 1 (and recall that we are still 
free to specify the value of M). In this case, we have 


2 1 


< 


2 1 


< 


21 


t 2 + u 2 ~ t 2 + M 2 t ~ M 2 t 


2 

M 2 ' 


Using (100), we see that if we assume M to be large enough, this forces 
y = 2t/(t 2 + u 2 ) to be such that J^{iy) < 1.01/yd/. Selecting such a value of 


M, and restricting /3 to be small enough to ensure that the (l + O (e 2,r /U) 


term in (101) is also < 1.01, we get therefore that the right-hand side of 


(101) is bounded by 


1.03 


t 2 \ 1/4 (f 2 + u 2 ) 1 / 2 _ 1.03 


t 2 + 


u- 


(2t)V2 


V2 


(t 2 + u 2 ) 1/A < 0.8(1 + t 2 ) 1/4 . 


If we assume further that /3 < 1/2, this is bounded by 0.9, which proves the 
claim in this case. 


Cases 1, 2 and 3 above cover all the possibilities. We have therefore 
proved Proposition! 10.2 which as we have seen implies Theorem 10.1 □ 
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Proof of Theorem 9.2. Decompose the integral in (87) as 


pDn 1 / 5 

(t>t n {u/a n )e~ inu,(yn du= / (ft n (u/a n )e~ mu/(Jn du 

J-Dn !/5 

+ [ ft n (u/a n )e~ inul ^ du. ( 102 ) 

J [—7rcr n ,7rcr n ]\[— Dn 1 / 5 ^n 1 / 5 ] 


Here, the first integral converges to ff° e “ 2 / 2 du = \/2n by Theorem 


9.3 


Lemma 9.6 and the dominated convergence theorem. The second integra. 


can be seen using Theorem 10.1 to satisfy the bound 


' [—7rcr n ,7rcr n ]\[— Dn 1 / 5 ^Dn 1 / 5 ] 


4>t n (u/a n )e mu/an du 


< 2na n exp 


(—£n 3 / 10 ) , 


and therefore converges to 0. It follows that the left-hand side of (102) 


converges to \/2tt. which is precisely (87). As we have seen, that result 
implies Theorem |9.2| □ 


11 Proof of Theorem 11.11 


Rewrite (74) in the form 

r(n) = G(e~ tn ) exp(nf n )P tll (!V = n). 


(103) 


The last factor on the right-hand side is given asymptotically by (84), so the 


proof of Theorem 1.1 reduces to writing down explicitly the asymptotics of 


the first two factors, using Theorem 6.2 and the definition of t n in (82), and 


then combining the pieces. First, we have that 

exp (nt n ) = exp (Vin 2 / 5 — T 2 n 3//1 ° — Tsn 1 ^ 5 ^ . 


(104) 


Second, to expand G(e~ tn ), recall that t n can be expressed in the form 
tn = R 9 6 (l — + vq)) as in (79), in terms of the parameters u = T 2 /T 1 , 

v = 73 / 71 , q = nT 1 / 10 defined in (78). We then have that 


G(e tn ) = exp ( mt n 2/3 + fJ> 2 t n L/ ' z - 5 log t n 


f-1/2 _ 1 


+ u/(0) + | log 2 + 0(re -3 ^ 10 )^ , (105) 
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and each of the non-constant summands in the expansion can be evaluated 
asymptotically using a Taylor expansion, similar to the derivation in Sec¬ 
tion [8| The computation is a bit more tedious, since this time we need to 
consider the expansion up to the fourth order to attain the desired precision 
up to a o(l) error. 

_ o /o - 

Start with the term proportional to t n . We have (compare with (80)) 

t- 2 / 3 ={r iq 6 )- 2/3 (l- q (u + vq)r 2/3 

= rr 2 / V 4 + ^ q(u + vq ) + ^- q 2 {u + vq ) 2 

1 2-5-8 3/ , 3 1 2-5-8 - 11 4 . .4 5 A 

+ 6 ‘3^ {U + Vq) + 24 ' 3-3-3-3 g {u + V<l) + ° ((/ ' } 


= T- 


-2/3 / -4 


g" 1 + ( 0“ ) q~ 3 + l n v + ) q 


21 -2 


10 2-5-8 , 

+ —uv -\ - u q 

9 6 • 3 • 3 • 3 1 y 


/5 2 2-5-8 2 2-5-8-11 4 » n _ , . 

+ l 9" + 6^TIT3 3 “° + 24.3.3.3.3 " )« +0(ll) f 


— 1/2 

Similarly, for the t n ' term we have 

in 1/2 = (^ 6 ) _ 1 / 2 ( i - g(u + ^))- 1/2 

= r[ l/2 q~ 3 ^1 + ^q(u + vq) + ^ • ^g 2 ^ + vg) 5 


1 1-3-5 


+ x • 


= T - 


—!/2 / -3 


g °+ I I g 2 +|-w + 


6 2 - 2-2 
1 1-3 


q 3 {u + vqf + 0{q A ) 


2 - 2-2 


^ 2 g- 1 


+ 


1-3 1-3-5 n _ . . 

-2uv + „ n „ u 3 g° + 0(g) . 


2-2-2 6 • 2 • 2 • 2 
The contribution of the logt n term is simpler to understand, since after 
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exponentiating it becomes 


exp (-^ logt n ) = 1/3 = {nQ 6 ) 1/3 (! + q(u + vq)) 1/3 


_ 1/3 —2 


= T- 


q-\l + 0(q)), 


and the error term 1 + 0{q) = 1 + 0{n 1 / 10 ), being a multiplicative term 
outside the exponential sum, is sufficient for our purposes. 


Substituting the expansions above into (105) and collecting terms in 
powers of n 1//10 , we conclude that G(e ~ tn ) has an asymptotic expansion of 
the form 

G{e~ tn ) = rr'/'V/ 5 exp (b^ + B 2 n 3 / 10 + B,n 1 ^ + T? 4 n 1/10 

+ B 5 + u/(0) + § log 2 + 0(n- 1/10 )) , (106) 

where Bi, B 2 , B%, i? 4 are constants, having the elaborate definitions 
e> - 2/3 

B 1 = /XlTj , 

E> 2 —2/3 . -1/2 

#2 = 3^ lT l “ + i i2T l • 

R 2 —2/3 . 5 —2/3 2 . 1 -1/2 

-63 = -Ml^i ^ + gMl^ U + 2^ T 1 

10 _2/3 40 — 2/3 3 , 1 - 1/2 . 3 — 1/2 2 

S 4 = y /ilTi UV + —M1T1 U + ~ll 2 r 4 V + -/X 2 T! U , 


R 0 -2/3 2 , 40 —2/3 2 . iiU -2/3 4 

B 5 = gPlT! V + —/iiTi « V + ^Vl T l u 

, 3 -1/2 5 —1/2 3 

+ -/X2T4 UV + — M2Tf U ' 

The definitions of B±, B 2 , B%, B^, B$ can be further simplified. The details 
of this purely routine algebraic simplification can be found in the companion 
package [03]. This results in the more pleasing expressions 

B\ = 3X 2 , 

B 2 = — -X^Y, 

10 

B 3 = - ijrV, 

3 100 


no 


li 


b 4 = — —x~ 7 y 3 , 

3200 

b 5 = —-—x~ 10 y 4 . 

2560 
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We are finally ready to derive ([3]). Multiplying (84), (104), and (106), we 
get using (103) that 


r(n) = (1 + o(l)) ^r, 1/13 ~^=X exp ( B 5 + u/(0) + | log2)^ n 3 / 5 

x exp ((Bi + Ti)n 2 ^ + ( B 2 - r 2 )n 3,/10 + ( B 3 - r 3 )n 1 / 5 + H 471 1 / 10 ) . 

This expansion is precisely of the form ([3]), with (one can easily check) our 
claimed values ([6])-([9]) for the constants A\, A 2 , A 3 , A 4 and the value (19) 
for I\, which in conjunction with Theorem 1.5 implies (10). The proof of 
Theorem 1.1 is complete. 0 


12 Final comments 


12.1 Additional summation identities 


Our proof of Theorem |1.6| relied on computer experimentation to discover 
the coefficients a n j in ( |30| ). (Moreover, as we commented in Section [2j the 
coefficients (3 nt k could also be guessed in principle although we happened 
to know them in advance.) This suggests the possibility of discovering and 
proving additional identities, and perhaps even automating the process lead¬ 
ing to such results. For example, by slightly generalizing the system of linear 


equations (29) to cover different mod 6 congruence classes, we were able to 


discover empirically the following additional identities: 


BQn-2 
6 n — 2 


1 f 471 ') V ( 2n ^ p(n k) B2n+2k ~ 2 B in—2k 

2(6 n + 1) \2nJ ^ \2k - l) [ ’ ’ 2n + 2k - 2 An - 2k' 

(107) 


B t 


6 n 


(4n + 1)! A ( 2 n 


6 n 


3(6n + 1) (2n)! 2 f^\2 k 




B 2n+2k B^ n _2k 

' 2n + 2k An — 2k' 


(108) 


where 


P(n, k ) = (2 n — l) 2 — A(k — l)(n — k), 

(re 2 — kn + fc 2 )(4(n 2 — kn + k 2 ) — 1) — 6 kn(n — k ) 


Q(n, k) = 


n( 2 k + l)(2n — 2 k + 1) 


Since the discovery was based on the summation technique of Sections [2] 
and[3j these identities should arise as specializations of similar identities for 
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the Eisenstein series. We did not work out the details of a proof, but our 
methods ought to apply. 


12.2 Open problems 

This work saw the fruitful meeting of ideas from complex analysis, number 
theory, asymptotic analysis, combinatorics, probability theory, and repre¬ 
sentation theory. Our results suggest several open problems that remain 
ripe for further investigation. Among them are the following: 


1. Find additional lower order terms in the asymptotic expansion for 
r(n). 

2. Prove analogous versions of our results for other connected simple 
compact Lie groups and simple Lie algebras. 


3. 


Find a direct proof of the result © that cj(— n) = 0 for integer n > 1 
without relying on the Bernoulli number summation identity (20), 


thereby deducing (20) in a new way. 


4. Use ideas from the theory of modular forms to reprove (23 ^ and to 
discover and classify other identities with a similar flavor. Do modular 
forms give additional insight into the properties and significance of the 
SU(3 ) zeta function cu(s)? 
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